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Research  Summary; 

Plasma-dynamic  thrusters  designed  for  high  speed  low  thrust  space  missions,  are 
nonlinear  distributed  parameter  systems  with  very  complex  dynamics.  The  main 
objective  of  this  study  was  to  analyze  MPD  thrusters’  controllability  and  observability 
without  the  need  to  solve  the  engines  governing  partial  differential  equations.  This 
analysis  falls  into  the  category  of  stabilizability  of  distributed  parameter  systems  (DPS), 
which  is  the  most  challenging  problem  for  the  control  community.  During  the  three- 
year  period  of  this  investigation,  several  new  concepts  were  derived  and  developed. 
Following  is  the  highlights  of  the  main  accomplishments  in  stabilizability  and 
observability  of  the  MPD  thursters. 

1.  Integration  of  advancements  in  the  mathematical  theory  of  partial  differential 
equations,  dynamical  systems,  and  advanced  feedback  control  theorems. 

2.  Derivation  and  application  of  Lyapunov  theorems  for  the  case  of  distributed 
parameter  systems. 

3.  Application  of  Lyapunov’s  stability  theorems  to  MPD  systems. 

4.  Derivation  of  the  characteristics  of  a  general,  nonlinear  hyperbolic  systems 
represented  by  partial  differential  equations,  and  their  relations  to  controllability  of 
the  MPD  systems. 

5.  Stability  analysis  and  stabilizability  of  MPD  thrusters. 

6.  Controllability  and  observability  of  distributed  parameter  systems. 

This  research  has  a  general  focus  on  the  analysis  of  the  DPS,  but  the  results  at 
every  stage  is  implemented  on  an  MPD  thruster  model.  Therefore,  the  new  foundations 
and  theorems  resulted  from  this  study  have  potentials  for  applications  in  many  other 
distributed  systems  than  MPD,  such  as  flexible  structures,  jet  engines,  integration  of 
controls  for  frame  vibration  and  combustion  process  in  an  aircraft  engine,  etc. 
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CHAPTER  1  -  INTRODUCTION 

Distributed  parameter  systems  (DPS),  also  regarded  as  infinite  dimensional 
systems,  are  often 'described  by  a  set  of  linear  or  nonlinear  partial  differential,  integro- 
differential  or  differential-delay  equations.  Conversely,  lumped  parameter  or  finite 
dimensional  systems  are  represented  by  a  set  of  ordinary  differential  equations.  For 
distributed  systems,  state  variables  are  functions  of  time  and  another  set  of  parameters, 
whereas  lumped  parameter  systems  are  only  dependent  on  time.  In  many  physical 
systems  the  states  of  DPS  are  functions  of  time  and  spatial  coordinates. 

In  general,  all  physical  systems  are  intrinsically  distributed  in  nature.  However,  in 
many  instances  the  system’s  spatial  energy  distribution  is  sufficiently  concentrated  such 
that  an  approximate  lumped  parameter  description  may  be  an  "adequate"  representation 
for  the  system.  On  the  other  hand,  the  spatial  energy  distributions  of  many  practical 
systems  are  widely  dispersed  and  require  controls  of  some  or  all  of  the  spatially 
distributed  states.  In  general,  such  systems  have  to  be  represented  by  distributed 
parameter  models.  Typical  examples  of  DPS  are  combustion  processes  in  engines, 
boilers,  and  furnaces;  heat  exchangers;  distillation  processes;  nuclear  and  chemical 
reactors;  gas  dynamic  thrusters  and  lightly  damped  flexible  structures.  Another  example 
of  DPS  with  a  widely  dispersed  distribution  of  electro  magnetic  energy  is  a  magneto¬ 
plasma  dynamic  accelerator.  Investigation  of  the  dynamics  of  distributed  parameter 
models  pertinent  to  magneto-plasma  dynamic  accelerators  is  of  interest  to  the  Air  Force 
for  future  space  vehicles  and  has  been  one  of  the  motivations  for  this  research. 

Often  a  nominal  (approximate)  model  for  a  system  is  considered  in  order  to 
provide  numerical  estimates  about  the  characteristics  of  DPS,  as  in  determining  the 
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system  response  to  an  excitation  and/or  design  a  control  input  to  achieve  a  desired 
system  performance.  These  approximate  models  are  often  derived  by  means  of  some 
modal  reduction  techniques.  The  most  common  methods  for  derivation  of  finite  order 
models  are  discretization  and  modal  truncation.  Due  to  the  reduction  of  the  system 
order,  some  of  the  information  about  the  dynamic  characteristics  of  the  system  is 
neglected.  This  implies  that  there  will  be  an  error  signal  between  the  exact  and 
approximate  solutions,  based  on  the  unmodeled  dynamics  of  the  system.  Usually,  if  an 
upper  bound  for  such  approximations  exists,  then  that  bound  grows  with  the  reduction 
of  system  order.  The  higher  the  approximation  bound,  the  lower  the  precision  in  the 
system  representation.  To  increase  accuracy,  the  number  of  discretized  points  (nodes) 
needed  may  reach  several  hundred.  Application  of  these  large  order  models  in  control 
design  would  require  a  large  number  of  sensors  and  actuators,  otherwise  interferences  ol 
uncontrolled  modes  would  affect  the  system  performance.  This  condition  is  often 
impractical  for  economical  reasons. 

One  of  the  main  concerns  of  control  theory  is  whether  or  not  a  system  is  stable. 
For  finite  dimensional  systems,  Lyapunov  stability  theory  have  become  an  important 
vehicle  in  the. stability  analysis  of  a  system.  This  approach  attempts  to  make  statements 
about  a  dynamic  system’s  stability  of  motion  without  explicit  knowledge  of  the  solutions 
to  its  governing  equations.  Although  the  development  of  Lyapunov’s  stability  theory  for 
ordinary  differential  equations  has  been  widely  investigated,  its  application  to  solutions 
of  partial  differential  equations,  namely  distributed  parameter  systems,  has  been  limited. 
Stability  results  for  distributed  parameter  systems  have  been  derived  based  on  spectral 
analysis  and  frequency  domain  characteristics  of  the  system.  Such  analysis  provides 
necessary  conditions  for  stability,  in  contrast  to  the  Lyapunov  method  that  results  in 
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sufficient  conditions.  In  order  to  obtain  sufficient  conditions  for  stability  based  on 
spectral  analysis  of  DPS,  other  system  properties  are  required.  Moreover,  frequency 
domain  analysis  of  DPS  is  often  used  to  provide  an  input-output  type  of  stability  which 
do  not  necessarily  imply  internal  stability  of  DPS.  These  requirements  are  in  addition  to 
the  complexities  associated  with  the  spectral  analysis.  In  general,  spectral  analysis 
requires  an  explicit  knowledge  of  solutions  to  system  equations.  Furthermore, 
approximate  models  cannot  be  used  to  imply  sufficient  conditions  for  stability  of  DPS. 
Considering  these  issues  indicates  that  the  Lyapunov  method  may  be  a  more  suitable 
alternative  for  stability  analysis  of  DPS. 

In  cases  where  system  performance  is  not  satisfactory  one  needs  to  determine 
whether  the  application  of  an  external  control  input  will  provide  the  desired  results.  To 
obtain  such  information,  the  system  must  be  "observable",  namely,  one  must  be  able  to 
estimate  its  characteristics  from  measurements  of  its  states.  Moreover,  the  system  must 
be  "controllable".  In  qualitative  terms,  based  on  the  knowledge  obtained  from 
observations,  one  must  be  able  to  manipulate  the  system  to  accomplish  a  desired 
performance.  These  concepts  for  finite  dimensional  linear  time  invariant  systems  have 
been  derived  and  are  well  known.  Extension  of  these  results  to  time  variant  and 
nonlinear  finite  dimensional  systems  can  be  established  in  the  local  sense  by  use  of 
implicit  function  theorem.  However,  the  technical  complexity  in  dealing  with 
controllability  and  observability  of  DPS  arises  from  the  fundamental  differences  between 
properties  of  finite  and  infinite  dimensional  spaces.  Hence,  the  approximate  models  have 
subtle  shortcomings  to  answer  questions  about  controllability  and  observability  of  DPS. 

In  cases  where  the  system  is  not  stable  or  asymptotically  stable,  one  must  know 
how  to  make  the  system  stable  or  asymptotically  stable.  This  process,  called 
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stabilization,  in  turn  depends  on  the  system  characteristics,  i.e.,  whether  or  not  it  can  be 
"stabilized/'  Although  this  property  for  DPS  is  related  to  observability/controllability, 
the  relationship  is  not  similar  to  the  way  these  three  properties  are  coupled  in  the  finite 
dimensional  systems.  Therefore,  the  approximate  and  reduced  order  models  of  DPS  are 
usually  insufficient  for  stability  analysis  and  the  design  of  stabilizing  controllers. 
Control  actions  based  on  such  models  may  have  interactions  with  uncontrolled  modes  of 
the  system  which  may  have  adverse  effects  on  stabilization  of  the  system. 

The  phenomenon  of  excitation  of  unmodeled  dynamics,  or  "spillover",  is  an 
example  of  such  adverse  effects.  Problems  of  this  nature  can  be  avoided  either  by  means 
of  distributed,  i.e.,  "body  force"  type  control,  or  by  finite  order  controls  based  on  a 
distributed  model  of  the  system.  In  applications  of  finite  order  (not  spatially 
distributed)  controls  where  actuators  are  concentrated  at  certain  finite  number  of 
locations,  stabilization  of  the  system  is  not  always  guaranteed  unless  the  points  of 
measurement  and  actuation  of  DPS  are  strategic  points.  Often,  the  boundary  of  DPS, 
or  parts  of  it,  has  strategic  point  properties.  This  leads  to  boundary  control  for 
stabilization  of  DPS.  These  issues  are  the  current  research  topics  in  the  control  area. 
This  study  attempts  to  develop  answers  for  these  questions  for  the  case  of  hyperbolic 
and  parabolic  types.  The  following  section  describes  some  of  the  previous  work  in  this 
area. 

Literature  Review 

The  first  significant  work  towards  control  of  DPS  was  initiated  by  Butkovskii  and 
Lerner  [1]*  Their  work  has  been  concentrated  on  the  derivation  of  maximum  principle 


-  8  - 


applied  to  a  class  of  DPS.  Subsequently,  [2|  studied  the  optimal  control  problem  of  a 
specific  linear  diffusion  system,  A  complete  description  of  all  relevant  issues  involved  in 
modeling,  stability,  control  and  optimization  of  DPS  was  first  reported  by  Wang  ;3i. 
The  studies  of  distributed  parameter  systems,  as  systems  of  partial  differential  equations, 
have  been  the  result  of  the  developments  in  the  mathematical  theory  of  partial 
differential  equations  (PDE)  and  dynamical  systems.  In  fact,  the  progress  made  in  the 
developments  related  to  DPS  are  due  to  the  concurrence  of  three  fields  of  science: 
controls,  PDE  and  dynamical  systems.  The  theory  of  PDE,  which  was  developed  to  deal 
with  problems  and  phenomenon  of  continua,  is  associated  with  the  names  of  many  great 
mathematicians.  More  modern  and  abstract  approaches  to  this  field  are  due  to 
Hadamard,  Lax,  Sobolev  and  Hilbert  [4-9].  The  theory  of  dynamical  systems  was 
initiated  by  the  pioneering  work  of  Zubov  [10|.  He  tried  to  make  mathematical 
abstraction  of  physical  systems,  called  dynamical  systems,  in  order  to  distinguish  general 
properties  of  these  systems,  such  as  stability  and  asymptotic  behavior.  With  this 
categorization  he  could  extend  stability  results  of  Lyapunov  [11]  to  systems  of  PDE. 
Sufficient  conditions  for  the  stability  of  equilibrium  solutions  for  a  system  of  PDE  were 
derived  by  Massera  [12].  The  application  of  Lyapunov  stability  theorem  based  on  the 
work  of  Zubov  has  been  investigated  by  several  authors.  Hsu  ]l3]  applied  this  theory  to 
a  nuclear  reactor  sy-tem.  Wang  [3]  considered  the  stability  of  those  evolution  equations 
whose  solutions  involve  a  semigroup  property.  There  are  also  many  other  applications 
which  utilize  Lyapunov  functions  directly  to  study  special  problems  [14,15,16^  .\ 

completely  rigorous  and  abstract  approach  to  the  theory  of  Lyapunov  stability  for 
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infinite  dimensional  systems  was  studied  by  [17]  and  [18].  Stability  study  of  an 
equilibrium'  solution  of  a  magneto-plasma  dynamic  (MPD)  system  for  the  special  case, 
where  the  plasms  equilibrium  velocity  is  zero  was  addressed  in  [19].  In  cases  where  a 
system  does  not  have  an  equilibrium  state,  the  invariance  principle  provides  information 
about  its  asymptotic  behavior.  The  invariance  principle  was  first  introduced  by  LaSalle 
[20]  for  finite  dimensional  systems.  The  generalization  of  this  principle  for  abstract 
dynamical  systems  has  been  accomplished  by  Hale  [21].  Applications  and  contributions 
in  the  asymptotic  behavior  of  an  abstract  d>..,amical  system  have  been  investigated  by 
[22-26]. 

Other  developments  in  the  theory  of  DPS  have  come  from  advances  in  the  control 
theory  applied  to  DPS.  Historically,  for  finite  dimensional  systems,  the  concept  of 
controllability  was  introduced  in  the  early  1960s  by  Kalman  [27,28],  The  concept  of 
controllability  was  extended  to  DPS  initially  by  Fattorini  [20-31],  who  investigated 
controllability  analysis  of  a  heat  equation.  Subsequently,  Triggiani  [35]  has  developed 
controllability  and  observability  concepts  for  general  DPS.  Applications  of  these  results 
to  the  hyperbolic  boundary  value  problems  are  reported  in  [36-40].  Other  applications 
of  observability  analysis  for  DPS  of  parabolic  type  are  reported  by  [41,43].  In  the 
process  of  the  evolution  of  controllability  concepts  for  DPS,  the  theory  of  time  optimal 
control  was  extended  to  DPS.  This  theory  was  developed  based  on  investigation  of 
many  researchers,  especially  the  work  of  Butkovskii  and  Lions  [44,45].  Lions’  work  is  a 
landmark  in  this  category.  Also,  with  more  emphasis  being  placed  on  hyperbolic 
boundary  value  problems,  investigations  of  [46,47]  can  be  recalled.  Moreover,  due  to  the 
nature  of  DPS,  the  optimal  control  problem  can  be  defined  in  terms  of  optimal  locations 
for  sensors  and  actuators.  Measurements  at  certain  points  in  the  spatial  domain  of  the 
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system  may  yield  more  information  about  the  system  than  other  points.  Since  the 
number  of  sensors  is  generally  governed  by  economical  considerations,  it  is  desirable  to 
locate  the  given  nuinber  of  measurement  sensors  at  points  that  lead  to  the  best  estimates 
of  the  system.  The  existence  theorem  for  the  solution  of  this  optimal  location  problem 
was  proven  by  Bensoussan  [48].  An  algorithm  for  derivation  of  suboptimal  sensor 
locations  based  on  eigenfunction  expansion  representation  of  linear  DPS  was  developed 
by  Yu  and  Seinfeld  [49].  The  dual  problem  of  optimal  actuator  location  for  DPS  has 
been  investigated  by  Amouroux  and  Barbary  [50]. 

Often,  a  control  problem  is  concerned  with  stabilization  as  well  as  the 
aforementioned  performance  optimization.  Research  on  the  problem  of  stabilizability  in 
infinite  dimensional  spaces  was  first  initiated  by  Slemrod  [51,52].  His  work  was 
primarily  motivated  by  hyperbolic  systems,  and  he  used  a  generalization  of  the 
invariance  principle  to  study  stabilization  of  infinite  dimensional  systems  in  Hilbert 
spaces.  This  problem  was  also  treated  using  the  spectral  theory  of  unbounded  operators 
for  stabilization  of  systems  in  Banach  spaces  by  Triggiani  [53].  Stabilization  of  DPS 
based  on  [52,53]  has  been  applied  to  several  classes  of  systems  [54-56]. 

Finite  order  control  problems  for  the  stabilization  of  DPS  has  been  studied  by 
Balas  and  Slemrod  [57-59].  Balas  [57]  studied  a  flexible  structure  with  pointwise  sensor 
and  actuator  to  control  a  modal  truncated  model  of  the  system,  based  on  the  original 
DPS.  The  truncated  model  is  not  always  controllable  and/or  observable.  It  is  only 
controllable  when  an  estimate  is  determined  to  prevent  spillover  problems.  This 
estimate  is  related  to  the  initial  perturbations.  In  [58,50]  boundary  control  of  the  wave 
equation  and  a  flexible  beam  subject  to  the  boundary  controls  were  investigated, 
respectively. 
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This  study  considers  stability  of  two  general  models  of  DPS  based  on  Lyapunov's 
direct  methods.  The  controllability  and  observability  concepts  based  on  abstract 
dynamical  system  are  analyzed  ana  their  relations  to  stabilization  by  distributed  and 
finite  order  controllers  are  derived.  Although  the  basics  of  stabilization  by  spectral 
theory  of  unbounded  operator  is  described,  the  fundamental  approach  for  stabilization  is 
based  on  [52]  and  the  invariance  principle. 
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CHAPTER  2  -  MATHEMATICAL  PRELIMINARIES 

This  chapter  provides  the  mathematical  bases  and  preliminaries  for  the  following 
chapters  of  this  study  to  provide  the  reader  with  a  quick  review  of  the  required 
definitions  and  theorems.  Additional  details  on  these  topics  can  be  found  in  references 
[60-65]. 

Vector  Spaces 

Definition  2.1.  A  real  (or  complex)  vector  space  is  a  set  X  such  that: 

1.  given  any  x,yeX,  there  is  an  element  x  -h  y  in  X  satisfying 

(a)  X  +  y  =  y  +  X ,  VA.yfX 

(b)  X  +  (y+z)  =  (x+y)  +  z ,  V/ic,y,zeX 

(c)  there  is  an  element  OeX  such  that  0  +  x  =  x ,  VAeX 

(d)  given  an  xeX,  there  is  an  element  — x  in  X  such  that  x  -|-  (— x)  =  0; 

2.  given  any  xeX  and  any  number  aeF  (F  being  a  real  or  complex  field),  there  is  an 
element  axeX  such  that 

(e)  a{i3\)  =  (o;/?)x,  for  any  a,f3eF  and  xeX 

(f)  (a-f.^x  =  ax  +  3x,  \/a,:3€F  and  xeX 

(g)  a(x+y)  =  ax  -h  ay,  V/afF  and  x,yeX 


(h)  lx  =  X,  VAeX 
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Definition  2.2.  A  subset  W  of  vector  space  X  is  a  subspace  if  it  is  a  vector  space  with  the 
operations  of  addition  and  multiplication  defined  on  X. 

Definition  2.3.  A  se't  of  vectors  {xi,  X2,  x^}  are  linearly  dependent  if  there  exists  a 
set  of  scalars  aj,  OfnCF  not  all  zero,  such  that  ajXi  +  <^2X2  +  •  •  •  +  =  0. 

If  no  such  set  of  scalars  exists,  then  Xi,  X2,  •••,  Xq  are  linearly  independent. 

Definition  2.4.  A  linear  combination  of  a  finite  set  of  vectors  Xi,  X2,  ...,  is  a  vector  of 
the  form  ctiXi  -t-  0:2X2  +  •  •  *  +  o^Xq,  where  Oj’s  e  F. 

Definition  2.5.  Let  W  be  any  set  of  vectors  in  a  vector  space  X.  The  set  of  all  linear 
combinations  of  elements  of  W  is  called  the  subspace  spanned  (generated)  by  W, 

Definition  2.6.  A  set  of  {xi,  X2,  ...,  Xq}  is  called  a  finite  basis  for  a  vector  space  X  if  it  is 
linearly  independent  and  it  spans  X. 

Definition  2.7.  A  vector  space  is  said  to  be  n-dimensional  if  it  has  a  finite  basis 
consisting  of  n  elements.  A  vector  space  with  no  finite  basis  is  said  to  be  infinite 
dimensional.  For  example,  the  space  of  continuous  functions  on  the  domain  ia,bj  (called 
C[a,b])  is  infinite  dimensional. 

Normed  Linear  Spaces 

Definition  2.8.  Let  X  be  a  vector  space  over  a  real  or  complex  field  F.  A  norm  on  X, 
denoted  by  II'  II,  is  a  real-valued  function  on  X  with  the  following  properties: 

a)  llx  II  >  0  if  X  ^0  and  llx  II  =  0  for  x  =  0,  for  all  xeX 

b)  llax  II  =  vbavb  llx  II,  for  oeF 


-  14  - 


c)  llx  +  y  II  <  llx  II  +  lly  II,  for  x  and  yeX 
Norms  can  be  constructed  in  different  ways. 


If  X  =  {xi  vb  i=l,2...,n.} 

then  for  infinite  dimensional  space  X,  n  — ►  oc,  the  following  is  defined: 


norm  =  II*  H/’p  = 


E 

i  =  l 


1/P 


<  oc,  and  P  >  1 


V  xf 
i-l 


1/2 


€2  norm  ^  II*  11^^  = 

^oa  =  11'  ll^*  =  Sup  vbXi  vb 


<  00 


L2  norm  ^  Hf  Hl,  = 


f  vbf(t)vb^  dt 


1/2 


,  f(t)efunctional  space  X 


Lpj  Loo  norm  are  defined  correspondingly. 


Definition  2.9.  A  linear  vector  space  with  a  norm  II*  \  is  called  a  normed  linear  space 
and  is  denoted  by  (X,  II*  llj). 

For  example,  the  sets  of  real  (and  complex)  numbers  R“  (and  C^)  are  linear 
normed  spaces  with  norm: 


lx  II  = 


5]  vbXi 


1/2 


V/6c€R“  (or  xeC°,  respectively) 


The  previously  mentioned  vector  space  C[a,bj  of  continuous  functions  on  the  interval 
[a,b]  can  be  turned  into  a  normed  linear  space  by  defining  the  norm 

1/2 


Ilf  II  = 


f  vb  f(t)  vb^  dt 


V/feC[a,b] 


or 
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llf  II  =  Sup  (vbf{t)vb  :  a<t<b}  V/feC[a,bj 

where  each  norm  can  be  proven  to  satisfy  the  properties  (a),  (b),  (c)  in  definition  2.8. 

Lemma  2.1.  For  vectors  x  =  {xj,  x^}  and  y  =  {yi,  yn}, 

vbf]Xiyivb*<  E^bXjVb*  =  llx  11^  lly  11^ 

i-l  i'-l 

where  y  denotes  complex  conjugate  of  y.  This  inequality  is  called  Cauchy-Schwartz 
inequality. 

This  same  inequality  is  true  for  infinite  dimensional  spaces 

b  _  b  b 

vb  f  f(t)  g(t)  dtvb*  <  (f  vbf(t)vb*  dt)  (f  vbg(t)vb*  dt) 

i  Z  A 

<  llf(t)ll*  llg(t)lP 

Definition  2.10.  A  subset  S  of  a  normed  linear  space  X  is  called  bounded  if  there  is  a 
number  M  such  that  llx  II  <  M  for  all  xcS. 

Definition  2.11.  A  map  F:  X — ►Y  between  two  normed  linear  spaces  (X,  II*  11^)  and 
(Y,  II*  lly)  is  said  to  be  continuous  at  x^eX  if  for  a  given  e  >  0  if  there  exists  a  5  >  0 
(5  =  5(e,Xo))  s^ch  that 

IIF(x)-F(xo)ll,  <e 

whenever  llx  —  x<,  II  <  ^. 

Definition  2.12.  A  sequence  {xq}  in  a  normed  linear  space  (X,  H*  ll*)  converges  to  Xq  if 

lim  llxn  —  Xq  llx  =  0  . 

n— *oo 

Continuity  and  convergence  are  closely  related  and  a  map  F  :  X  — Y  is  continuous  if 
and  only  if  (iff) 
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lim  F(Xq)  =  F(  lim  Xq) 

n— *00  n—^oo 

Example  2.1.  The  sequence  {f“n(*')}.  where  =  e““*  and  fu  e  C[0,  ij,  converges  to 


0  for  t  ^  0 


F{t)  = 


1  for  t  =  0 

in  both  Loo  L2  norms.  It  is  clear  that  F(t)  ^  C[0, 1].  If  one  considers  Fi(t)  =  0  for 
all  t  [0,1],  then  {fn(t)}  converges  to  Fi(t)eC[0,  l]  in  the  L2  norm  but  {fn}  does  not 
converge  to  F 1  (t)  in  Loo  norm. 


Definition  2.13;  Open  and  Closed  Set.  A  set  A  in  a  normed  linear  space  X  is  closed  if 
all  convergent  sequences  in  A  have  their  limit  points  in  A. 

A  set  A  is  open  if  its  algebraic  complement  is  closed.  Alternatively,  a  set  A  is  open 
if  for  any  point  xeA,  there  is  e  >  0,  such  that  the  set  {y :  llx— y  II  <  e}  is  wholly  contained 
in  A.  The  closure  of  a  set  can  be  formed  by  adding  all  limit  points  of  sequences  in  A  to 
A,  and  is  denoted  by  A. 


Definition  2.14:  Cauchy  Sequence.  A  sequence  {xn}  of  elements  in  a  normed  linear 
space  (X,  II  •  ll,c)  is  termed  Cauchy  if 

llXn  —  Xni  llx  — ►  0  as  m,  n  — ►  cc  . 

Definition  2.15.  A  normed  linear  space  is  complete  if  every  Cauchy  sequence  defined  on 
that  space  converges  to  a  limit  point  in  the  space.  A  complete  normed  linear  space  is  a 
Banach  space. 

Theorem  2.1.  The  space  C[a,b]  with  norm  Ilf  I!  =  Sup  {vbf(t)vb  :  a  <  t  <  b}  is 
complete. 
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The  space  of  continuous  function  with  the  square  integral  norm  (L2)  is  an 
incomplete  space.  This  fact  is  clear  from  example  2.1.  The  space  of  L2[a,b)  (square 
integrable  functions)  provides  the  completion  of  C[a,bl  with  respect  to  L2  norm.  Thus, 
Ljfajbj  contains  all  functions  which  are  the  limits  of  continuous  functions  in  the  sense  of 
mean  square  (L*)  convergence.  This  property  of  L*[a,b]  will  be  more  elaborated  under 
the  topic  of  distributions  and  Sobolev  spaces. 

Definition  2.16.  A  subspace  S  of  a  normed  linear  space  X  is  called  dense  in  X  if  its 
closure  with  respect  to  the  norm  is  equal  to  X. 

Equivalent  Norms 

Definition  2.17.  Let  11*  llj  and  II  II2  be  two  different  norms  on  the  same  vector  space  X. 
A  norm  II  llj  is  called  equivalent  to  II*  II2  if  there  are  positive  numbers  a  and  b  such  that 

a  llx  111  <  11x2  II  <  b  llx  111  ,  \/ xeX 

It  is  easy  to  show  that  if  II*  lli  is  equivalent  to  II*  Il2,  then  II*  II2  is  equivalent  to  II*  llj 
with  bound  coefficients  replaced  by  b“^  and  a”^. 

Theorem  2.2.  If  a  sequence  converges  with  respect  to  one  norm,  then  it  converges  with 
respect  to  any  norm  equivalent  to  it. 

Corollary  2.1.  If  a  space  is  complete  with  respect  to  one  norm,  then  it  is  complete  with 
respect  to  any  norm  equivalent  to  it. 

Theorem  2.3.  In  a  finite  dimensional  space,  all  norms  are  equivalent. 

Operators  on  Vector  Spaces 
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Definition  2.18.  Let  X  and  Y  be  subsets  of  two  vector  spaces  M  and  N,  respectively.  An 
operator,  mapping  or  transformation  T  ;  X  — ►  Y  is  a  rule  which  associates  any  given  xeX 
with  an  element  of  Y  and  is  denoted  by  Tx  or  T(x).  This  mapping  may  also  be  implied 
by  X  vb— ►  Tx,  which  implies  "x  is  mapped  into  Tx". 

The  sets  X  and  Y  are  called  the  domain  and  range  of  T,  respectively.  The  vector 
Tx  is  called  the  image  of  x  under  T. 

Definition  2.19.  A  mapping  T:X— ♦X,  where  X  is  a  subset  of  a  normed  linear  space  N, 
is  called  a  contraction  mapping  if  there  is  a  positive  number  a  <  1  such  that 

II  Tx  —  Ty  11  <  Q:  llx  —  y  II  for  all  x,yeX 

Theorem  2.4.  If  T ;  X  — ►  X  is  a  contraction  mapping  of  a  closed  subset  X  of  a  Banach 
space,  then  there  is  exactly  one  xeX  such  that  Tx  =  x.  The  point  x  in  this 
transformation  is  called  the  fixed  point.  This  theorem  gives  proof  of  existence  and 
uniqueness  and  a  computational  method  for  finding  the  solution  (fixed  point)  of 
differential  and/or  integral  equations. 

Compactness 

A  bounded  region  in  a  finite  dimensional  space  has  a  finite  volume.  Successive 
applications  of  an  operator  T  to  a  point  provides  a  sequence  of  infinite  points.  Since 
these  infinite  number  of  points  are  contained  in  a  bounded  volume,  then  the  sequence  of 
operation  converges  to  a  limit  point,  which  is  the  fixed  point  of  T.  But  in  an  infinite 
dimensional  space,  even  a  bounded  region  is  so  large,  due  to  the  infinite  number  of 
dimensions  into  which  it  extends,  that  a  sequence  can  wander  through  the  region 
indefinitely.  Therefore,  the  sequence  never  converges  and  never  reaches  a  fixed  point.  In 
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infinite  dimensional  spaces  compactness  is  used  to  define  boundedness. 

Definition  2.20.  A  subset  S  of  a  normed  linear  space  X  is  compact  if  every  infinite 
sequence  of  elements  of  S  has  a  sub-sequence  which  converges  to  an  element  of  S. 

In  finite  dimensional  spaces,  closed  and  bounded  sets  are  the  same  as  compact  sets. 
However,  closed  and  boundedness  are  not  sufficient  to  give  compactness  in  the  infinite 
dimensional  spaces. 

Linear  Transformations 

Definition  2.21.  A  linear  transformation,  T,  from  a  linear  space  X  to  a  linear  space  Y 
over  F  is  a  map  T  :  X  — ►  Y,  such  that 

T(q!X  +  dy)  —  Q;Tx  +  /?Ty 

for  all  x,yeX. 

Definition  2.22.  T  :  X  — ►  Y,  were  X  and  Y  are  normed  linear  spaces,  is  said  to  be 
bounded  if 

IITxll<Kllxllx 

for  some  constant  K  >  0  and  all  xeX. 


Definition  2.23.  If  T  is  a  bounded  linear  transformation  between  two  normed  linear 
spaces  X  and  Y,  one  can  define  II T  II  by 


IIT  II  =  Sup 

x^O 

xfX 


IITxl^ 

ITT 


Definition  2.24.  If  X  and  Y  are  normed  linear  spaces,  then  one  can  define  the  space  of 
all  bounded  linear  transformations  T  :  X  — *>  Y  by  if  (X,Y)  with  the  norm  defined  in 
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Definition  2.23. 

Definition  2.25.  If  T  is  a  linear  operator  (transformation)  defined  on  its  domain 
D(T)  C  X  with  range  R(T)  C  Y,  then  the  graph  G{T)  is  the  set 


G(T)  =|(x,  Tx),  xeD(T),  TxeR(T)] 
in  the  product  space  X  X  Y.  The  linear  operator  T  is  closed  if  its  graph  G(T)  is  a  linear 
subspace  of  XXY. 


Definition  2.26.  A  linear  transformation  over  pace  X,  which  maps  it  into  F,  is  called  a 
linear  functional 

T:X  — F 

Theorem  2.5  (Hahn-Banach  Theorem).  Every  continuous  linear  functional  h  :  M  — ♦  F 
defined  on  a  linear  subspace  M  of  a  normed  linear  space  X  can  be  extended  to  a 
continuous  linear  functional  H  on  all  of  X  with  preservation  of  norm. 

A  useful  corollary  of  this  theorem,  which  will  be  used  in  the  controllability  and 
observability  analysis  of  the  system,  is  given  as  follows: 


Corollary  2.2.  If  E  is  an  arbitrary  subset  of  a  normed  linear  space  X,  then  span  E  =  X  if 
and  only  if  the  zero  functional  is  the  only  functional  which  vanishes  on  all  of  E. 

Definition  2.27.  Dual  (conjugate)  space  of  a  normed  linear  space  X  is  the  normed  linear 
spree  of  all  bounded  linear  functionals  on  X.  One  can  write  this  space  as  X  (X,F)  or  X*. 

Definition  2.28:  Inner  Product.  An  inner  product  on  a  linear  space  H  defined  over  the 
complex  or  real  field  F  is  a  map  <-,*>:HXH— ►F  such  that 
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(a)  <Qac  +  /?y,  z>  =  a<x,z>  +  f3<y,z> 

(b)  <x;y>  =  <y,x> 

(c)  <x,x>  >  0  and  <x,y>  =  0  if  X  =  0 
for  X,  y,  z,  c  H,  and  a,  /3  eF. 

A  linear  space  H  with  an  inner  product  <•,*  >  is  called  an  inner  product  space. 

Definition  2.29.  A  Hilbert  space  is  an  inner  product  space,  which  is  complete  as  a 
normed  linear  space  with  the  norm  defined  by  inner  product. 

Ilx  IIh  =  <x,x> 

Definition  2.30.  An  orthonormal  space  of  a  separable  Hilbert  space,  H,  is  maximal  if 
H  =  span  {0i }.  Then  for  any  xeH,  the  Fourier  expansion  is 

X  =  S  <x. 

i»l 

One  of  the  most  important  properties  of  a  Hilbert  space  is  the  simple 
representation  of  its  dual.  The  following  theorems  illustrates  this  representation. 

Theorem  2.6;  (Riesz  Representation  Theorem).  Every  bounded  linear  functional  in  a 
Hilbert  space  H  can  be  written  in  the  form  <x,yo>  where  yoeH  is  uniquely  determined 
from  T(x). 

T(x)  =  <x,yo  >  ,  VAeH 

Theorem  2.7:  (Lax-Milgram  Theorem).  Let  H  be  a  Hilbert  space  and  let  B(x,y)  be  a 
complex  valued  functional  defined  on  the  product  Hilbert  space  HXH,  satisfying  the 
following  conditions: 
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i.  Sesqui-linearity,  i.e., 

B  [(q:iXi  +  02X2), yj  =  aiB(xi,y)  -h  a2B(x2,y) 

and 

B  [x,{/?iyi  +  52y2)j  =  A  B(x,  yO  -f-  ^2  B(xi,y2) 

where  0i,  02  complex  conjugates  of  0i,  82  respectively. 

ii.  Boundedness,  i.e.,  there  exists  a  positive  constant  7  such  that 

vb  B(x,y)  vb  <  7  llx  lb  Ny  II 

iii.  Positivity,  i.e.,  there  exists  a  positive  constant  S  such  that 

B(x,x)  >  (5 llx  iP 

Then  there  exists  a  uniquely  determined  bounded  linear  operator  Sei?{H,H)  with  a 
bounded  linear  inverse  S~*eif(H,H)  such  that 

<x,y>  =  B(x,Sy);  IISII<l/(5 

and 

<x,S“^y>  =  B(x,y);  IIS~Ml<7 

This  theorem  can  be  used  in  relation  with  equivalent  inner  products  and  norms. 

Theorem  2.8.  Two  inner  products  defined  on  a  real  linear  vector  space  H  are  equivalent 
if  and  only  if  there  exists  a  symmetric  bounded  positive  definite  linear  operator 
Sei?{H,H)  such  that 

<x,y>2  =  <:x,Sy>i 

where  the  indices  identify  the  inner  product  in  Hilbert  spaces  1  and  2.  Because 
B(x,y)  =  <x,Sy>i,  it  can  be  concluded  that  <x,y>2  is  sesqui-linear,  bounded  and 
positive  due  to  the  properties  of  S,  and  that  <x,y>2  satisfies  all  the  properties  of  an 
inner  product.  Hence  the  norms  defined  by  these  inner  products  are  equivalent.  The 
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application  of  equivalent  norms  yields  the  derivation  of  the  stability  properties  of  the 
general  operators  (A)  in  differential  equations  of  the  form 

X  =  Ax , 

from  the  knowledge  of  the  properties  of  the  operator  A,  The  concept  of  equivalent  inner 
products  can  be  extended  to  complex  Hilbert  spaces. 


Definition  2.31.  Let  A  e  i?  (H),  then  the  adjoint  operator  A*  is  defined  by 

<Ax,y>  =  <x,A*y>  for  all  x,  y  e  H  . 

This  suggests  that  for  bounded  operator  A,  the  adjoint  is  a  transpose  of  its  complex 
conjugate.  For  unbounded  operators,  where  x  e  D(A)  and  y  e  D(A*)  with  D(A)  being 
dense  in  H,  it  can  be  shown  that  if  A  is  closed,  then  D(A*)  is  dense  in  H  and  A*  is 
closed.  Hence 


For  example,  if 


A**  =  A 


Au  =  €  H  =  L2(0,T)  ,  with  u(0)  =  0 

T  . 

<Au,v>  =  f  v(t)  dt 

T 

=  u(T)  v{T)  -  /  u(t)  ^  dt 


then 


du 


where  u,  and  —  e  H  with  u(T)  =  0  . 

(I  w 


A*u  =  - 


du 

dT  ’ 


Distributions 

Definition  2.32.  The  support  of  a  function  f:R— *-C  is  the  set  {x:f(x)  #  0},  written  as 
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supp(f).  A  function  has  bounded  support  if  there  exists  numbers  a,b  such  that 
supp(f)  C  (a,b]. 

Definition  2.33.  Letting  0  be  an  open  set  in  R“,  then  the  function  f;!!  — ►  C  is  n  times 
continuously  differentialable  if  its  first  n  derivatives  exist  and  are  continuous.  The 
function  f  is  denoted  by  feC“(n).  f  is  smooth  or  infinitely  differentiable  if  f  e  C^.  The 
set  C5(fl)  consists  of  functions  in  €“(0),  which  vanish  outside  a  compact  subset  fl'  C  ^1. 

Definition  2.34.  A  test  function  is  a  smooth  R  — ►  C  function  with  bounded  support. 
The  set  of  all  test  functions  is  called  2. 

Definition  2.35.  A  linear  functional  on  2l  is  a  map  f :  ©  — ►  C  such  that 
f(a(^  +  htli)  =  a  f((^)  +  bf(t)  for  all  a,beC  and 

Definition  2.36.  A  sequence  of  test  functions  converges  to  a  limit  point  p  in  S,  if  (i) 
there  is  an  interval  [a,b]  containing  supp(p)  and  supp(0n)  for  all  n,  and  (ii)  for  any  k, 
0W(x)  0('‘)(x)  as  n  — ►  oo,  uniformly  for  x£[a,b],  where  is  the  k-th  derivative  of  d. 

Definition  2.37.  A  functional  f  on  is  continuous  if  it  maps  every  convergent 

sequence  in  2:{Q)  into  a  convergent  sequence  in  C;  that  is,  f(  Oq)  — *■  f(C')  whenever 
— *•  0  in  A  continuous  linear  functional  on  Q  is  called  a  distribution,  namely 

f  =  f(x)  o(x)  dx  =  <f,P>  =  f(p)  V/0€2(P.) 

The  set  of  all  distribution  in  the  compact  support  fl  is  denoted  by  P  (P). 

Definition  2.38.  Two  distributions  f  and  g  are  equal  on  (a,b)  if  <f,  o>  =  <g,  c*>>  for  all 
such  that  supp(<;^)  C  (a,  b). 
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The  delta  distribution,  5,  is  defined  as 

<S,  <f>>  =  (^o)  for  all  0  in  ^  . 

Definition  2.39.  The  derivative  of  a  distribution  f  is  a  distribution  defined  by 
<f^,0>  =  — <f,0*>  for  all  0e3. 


Similarly,  higher  order  derivatives  can  be  written  as 

<D“f,o>  =  (-1)”  <f,D“0>  . 

A  formal  partial  differential  operator  defined  in  an  open  subset  H  of  R“  can  be 
shown  in  abbreviation  by 


7"  =  E  aj  (x) 

vb  J  vb  <  m 

where  the  highest  order  of  differentiation  is  m,  symbol  J  denotes  an  index,  i,e.,  a  k-tuple, 

k 

J  =  (ji.  h>  -M  jk),  and  vbJvb  =  X)  Ji  , 

i-l 


= 


Jk 


The  adjoint  operator  of  r  is 


5x1*  5x2^  ’  •  ’  dx 


=  E  (-1)'^  ^  [aj(x) 

vb  J  vb  <  m  ^ 


Definition  2.40.  Let  7  be  a  formal  partial  differential  operator  defined  in  open  set  Q  e  R° 
with  real  coefficients  aj(x)  e  C^(n).  If  fe©'(n),  then  rf  will  denote  a  distribution  defined 
by  the  following  equation 

(rf)  (0)  =  f(r*<i))  ,  0€Cf(n) 


Sobolev  Spaces 


Definition  2.41.  Letting  H  be  an  open  subset  of  R“  and  let  K  be  a  non-negative  integer, 
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then 

(i)  the  set  of  all  distribution  feQ'  such  that  d^f6L2{^)  for  all  vbJvb  <K  will  be 

denoted  by  This  space  is  called  Sobolev  space,  which  is  an  inner  product 

space  mapping  each  pair  f,g 

<f.g>K  =  S  /  dx 

vbJvb  <  K  n 

The  norm  based  on  this  inner  product  is 

tIfllK  =  (<f.f>K)'^" 

(ii)  the  space  H^(n)  is  the  closure  of  C*^(f2)  in  the  above  norm,  similar  to  the  space 
H  *  (fl)  =  L2(n),  which  is  the  closure  of  C(n)  in  II  •  IIl2  norm. 

(iii)  the  symbol  H^(n)  will  denote  the  closure  of  (H)  in  II’  IIhk  norm. 

Lemma  2.2:  Letting  Q  be  an  open  set  in  R'^,  then  the  space  H^(n)  of  the  preceding 
definition  is  a  complete  Hilbert  space  and  the  space  Ho  is  a  closed  subspace  of  H^(n). 
In  addition, 

HS(n)  c  H°(n)  =  L^n) 

H^+^(n)  CH^(n)  K>o 

H^^‘(n)CH^{n)  K>o 

Theorem  2-9.  Letting  H  be  a  bounded  open  set  in  R“  and  dfl,  the  boundary  of  D,  be  a 
smooth  surface  with  no  interior  point  to  Q.{  =  Vt^dU),  then  H“(n)CC^(n),  where  m 
and  K  are  integers  with  m  >  K  +  n/2. 

Theorem  2.10.  If  fl  is  as  given  in  the  preceding  theorem,  then  the  mapping 
H“(n) — ►H®~^(n)  is  compact. 
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CHAPTER  3  -  DYNAMICAL  SYSTEMS  AND  LYAPUNOV  THEOREMS 
Distributed  Parameter  System  (DPS) 

In  this  section  DPS  and  lumped  parameter  systems  are  compared  and  general 
properties  of  these  systems  are  described.  The  dynamics  of  distributed  systems  are  often 
expressed  by  a  set  of  partial  differential,  integral  or  differential  delay  equations,  as 
opposed  to  lumped  parameter  systems  which  are  described  by  a  set  of  ordinary 
differential  equations.  For  distributed  systems,  the  state  variables  are  functions  of  time 
and  another  set  of  parameters,  whereas  lumped  parameter  systems  are  only  dependent 
on  time.  In  many  physical  systems  the  states  of  DPS  are  functions  of  time  and  spatial 
coordinates.  Fig.  3.1  illustrates  the  analogy  between  DPS  and  lumped  parameter 
systems.  In  this  figure,  at  each  instance  of  time,  the  state  of  a  lumped  parameter  system 
is  veR*.  However,  the  state  of  DPS  is  a  vector  function  of  one  spatial  coordinate,  x. 
Therefore,  each  of  the  two  states  of  the  DPS,  shown  in  Fig.  3.1,  consists  of  infinite  points 
and  belongs  to  an  infinite  dimensional  functional  space.  This  example  reveals  the 
infinite  dimensional  nature  of  the  DPS  as  opposed  to  finite  dimensionality  of  lumped 
parameter  systems. 

Semigroups  and  Abstract  Evolution  Equations 

The  simplest  finite  dimensional  systems  are  linear  autonomous  systems  which  can 
be  formulated  by 

v=Av  veR“  (3-1) 

v(to)  =  Vo 

where  v:  [o,T]  — ►R“  is  the  solution  of  the  system.  The  operator  (matrix)  A  is  bounded, 
Aeif(R®,R'‘).  Solution  of  the  homogeneous  system  (3-1)  is 


k  V.  ;t,OI  • 

V(g.x)  - 

\ 

.  V 

A 

V^(x) 

V 

Vjtt.o) 


•  V 

o 


v(t)  ■  ♦ 

o 


p. 


VjCt,/) 


vev  (Co^),  E^) 


T(t)  t  v^(x)_^  v(t,x) 
o 


v(t,X)  ■  T(t)  ♦  T(t-S)iU 
o 


Figure  3.1  Comparison  of  finite  and  infinite  dimensional  systems. 
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V  =  6'^^V 


which  satisfies 


=  Ae''-‘vo  -  Av 


where  e"^*  is  the  transition  matrix  defined  by, 


=  I  +  At  + 


aH^ 

2! 


+  ....  +  ^  + 

n! 


Usually,  external  effects  such  as  body  forces  can  be  included  in  the  above  system 
representation  by  a  vector  valued  function  g :  te[0,T]  — 


V  =  Av  +  g{t)  v(0)  =  Vo  (3-2) 

The  solution  is  then 

t 

v(t)  =  e'^^Vo  +  /  g(s)ds  (3-3) 

O 

The  excitation  term  g(t)  can  be  represented  in  terms  of  a  vector  valued  (input)  function 
u(t):  [0,T]  — with  g(t)  =  Bu(t)  where 

In  the  case  of  distributed  parameter  systems,  the  mathematical  description  is 
usually  given  by  a  set  of  partial  differential  equations.  In  order  to  generalize  the  results 
from  the  finite-dimensional  systems,  the  set  of  partial  differential  equations  can  be 
transformed  into  an  abstract  form  of  equations  (3-1)  and  (3-2). 

v=Av  veV  |[0,^] ,  E“j  (3-4) 

v(to)  =  Vo 

This  formulation  of  the  system  is  called  the  evolution  or  the  state  space  equation.  The 
state  space  of  the  system  V  belongs  to  some  Banach  space  or  more  commonly,  to  a 
Hilbert  space.  The  evolution  operator  A  in  equation  (3-4)  is  applied  to  the  v  and  it  may 
be  a  bounded  or  unbounded  operator.  An  example  of  bounded  operator  is  Fredholm 


-  30- 


integral  operator  for  which, 

Av  =  /  a(x,C)v(^,t)d^  (3-5) 

a 

where  a(x,^  is  the  kernel  of  the  operation.  The  unbounded  operator  is  often  in  the  form 
of  a  derivative  action  which  could  map  a  bounded  variable  into  an  unbounded  one.  For 
example,  the  wave  equation 


can  be  written  as 


at^ 


ax2 


+  f(x.t) 


(3-6) 


avi 

~di 


=  Vj 


(9v2  a*vi 


=  0* 


di  ax* 

Then  the  evolution  operator  A  would  be 


+  f(x,t) 


(3-7) 


A  = 


a* 

which  is  an  unbounded  operator  because  the  — —  is  unbounded. 

ax^ 

The  solution  of  the  evolution  equation  (3-4)  for  infinite  dimensional  systems  can  be 


given  as 


v(t,x)  =  T(t)  Vp  (3-8) 

where  the  operator  T(t)  (strongly  continuous  semigroup  or  group)  is  an  abstraction  of 

operator  e^‘  with  similar  properties  [66],  namely: 

(a)  v(0;Vo)  =  Vq  (i.e.,  at  time  0,  initial  state  is  v,,) 
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(b)  v(ti  +  ta;  Vo)  =  V  (ti;  v(t2,‘  Vo)]  =  v  (tj;  v{ti;  Vo)j 

(c)  v(t,Vo)  is  continuous  in  t  and  Vq  (i.e.,  if  one  of  them  changes  slightly,  then  the 
solution  should  not  be  changed  drastically). 

Based  on  abstraction  of  the  above  conditions,  the  following  definition  can  be  given: 

Definition  3.1.  A  Cq  (strongly  continuous)  semigroup  is  an  operator  T(t):  R  "  (V), 

where  .2*  (V)  is  the  space  of  ail  bounded  operators  on  V  into  V.  The  following  properties 
characterize  a  semigroup: 

(a)  T(0)=I 

(b)  T(ti  +t2)=T{ti)T(t2)  ti,t2>0 

(c)  lim  T(t)v  =  V,  for  all  veV  (i.e.,  T  is  strongly  continuous  at  t  =  0). 

The  operator  T(t)  will  be  defined  as  a  C<,  -  group  provided  t  e  (— oo,  oo). 

Fig.  3.1  depicts  the  similarity  between  e^^  and  T(t)  in  deriving  the  finite  and 
infinite  dimensional  systems  from  an  initial  state  v  to  a  final  state  v(t).  However,  the 
operator  T(t)  is  applied  to  the  state  v  instead  of  mere  multiplication.  Two  examples  of 
semigroups  are  given  here.  These  semigroups  are  generated  from  the  evolution  equations 
below. 

For  the  diffusion  process  in  heat  or  mass  transfer,  the  governing  dynamics  is  given 
by 

+  f(x,t)  ,  v(0,x)  =  Vo(x) 

The  evolution  operator  A  in  this  case  is 


(3-9) 
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A  =  a 


dx^ 


By  applying  the  Laplace  transform  to  the  homogeneous  evolution  equation,  one  can 
determine  the  semigroup  T(t)  generated  by  this  system.  If  v{s,x)  is  a  Laplace  transform 
of  v(t,x),  then 

2  d^v 
— *  _ 


sv  -  Vo  =  a* 


dx^ 


(3-10) 


sv  -a^-^  =Vo(x) 


if  v\^  is  the  green’s  function  for  the  system 


/  Vh(s,0  Vo(x-0  d^  =  /  vjj  ^3,(x-o]  Vo(0 


(3-11) 


svji  —  a* 


5®Vh 


ax* 

Vt  (H-oo)  =  Vii(-oo)  =  0 
.  vTvbxvb 

Vh  =  C 


lim  —a*  ■  ^  vb^lt  =  +a*  C 
el— 0  ax 


f  vj,  dx  =  0 

— VT  vbxvb 


VT 


vb(  -t- 


VT 


—  VT  vb  X  vb 


vb_. 


=  -}*  2a  C 


v7  ^ 

- e  ‘ 


=  +  2a^C^  =  l 
a 


C  = 


2a  "X/s 


Vh  = 


2a 


VT 


exp 


VT 


vbxvb 


(3-12) 


From  equation  (3-11),  v  can  be  calculated 
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1  r 

v  =  — _  /  exp 

I  V  S  -oo 


2a 


ir^ 


V7 


vbx— ^  vb 


vo(ad^ 


1  7 
rr-  / 


2a Vs  - 


v7 


vbx— ^vb 


vo(e)de 


V  =  /  d^Vo{^)ir^ 


2a 


VT 


exp 


VT 


vb  x— ^  vb 


The  inverse  Laplace  transform  is  [67|: 

/ 

Therefore,  for  X  =  — ^ 


exp 


-X74t) 


(3-13) 


(3-14) 


V  = 


1 


2a  V^ rrt  -oo 


/  exp 


(x-O^ 


4a*  t 


Vo(0 


Hence,  T(t)  will  be 


T(t)  [i  =  — -^7—  /  exp 
2aV7rt  -00 


4a*  t 


(i  de 


(3-15) 


(3-16) 


The  semigroup  for  the  first  order  wave  equation  can  be  derived  in  a  similar  way: 


dv  ^  ,  F 

A  <5(-) 

A  =  -u 

ax 


v(o,x)  =  Vo(x) 


(3-17) 


Similar  calculations  results  in  the  following  form  for  the  state  of  the  system  v  at 


time  t: 


v(t,x)  =  T{t)vo(x)  =  Vo(x-ut)  (3-18) 

Theorem  3.1.  Letting  T(t)  be  a  semigroup  on  to  X  (V),  where  V  is  a  Banach  Space, 


then 
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(a)  IIT(t)  II  is  bounded  on  every  compact  interval  of  l0,ocj  such  that  IIT(t)  II  <  Me*^  for 
some  M,  uje^. 

(b)  T(t)  is  strongly  continuous,  on  [0,ooj. 

The  proof  of  this  theorem  is  given  in  [66|. 


Definition  3.2.  The  infinitesimal  generator  of  a  strongly  continuous  semigroup  T(t)  is 
defined  by  A  when 

Av  =  lim  —  (T(t)v— v)  • 

t— 0"^  t 

V  } 

veD(A)CV 

It  should  be  noted  that  in  general,  A  will  be  an  unbounded  operator.  The  operator  A  is 
closed,  i.e.,  its  range  and  domain  converge  to  some  element  of  their  respective  spaces.  In 
addition,  the  closure  of  the  D(A)  covers  the  space  V,  namely  it  is  dense  in  V. 


Theorem  3.2.  Let  T(t)  be  a  strongly  continuous  semigroup  on  a  Banach  space  V  with 
infinitesimal  generator  A.  If  v,,  e  D(A),  then 


(a)  T(t)vo  eD(A)  for  t  >  0 

(b)  ■^(T(t)vo)  =  A(T(t)v<,)  =  T(t)Avo,  t  >  0 

W  -^(T(t)vJ  =A“(T(t)v,)  =T(t)A“v„  for  v,eD(A“);  t  >0 

t 

(d)  T(tK  -  Vo  =  -  /  T(s)Av„  ds;T  >  0 


(e)  D(A“)  is  dense  in  V  for  n  =  1,  2,  •  •  •  and  A  is  closed. 
The  proof  is  given  in  [62]. 
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For  finite  dimensional  systems,  the  Laplace  transform  of  is  L(e'^‘)  =  (SI— A)  This 

can  be  generalized  to  semigroups  by  the  following  proposition: 

Proposition  3.1.  If  T(t)  is  a  strongly  continuous  semigroup  with  infinitesimal  generator 
A,  then  ReCs)  >  ui  for  sep(A).  Where,  p(A)  is  defined  as  the  resolvent  set  of  A,  namely 


p(A)  = 


S:(SI-A)-‘  fif  (V);  bounded  linear  operator  on  V 


Note  that  w  is  given  by 


=  Inf 


cj: 


:  tlT(t)  II  <  Me-'S-  M, 


where  Inf  is  greatest  lower  bound  of  the  set  of  a;’s. 


Theorem  3.3:  (Hille-Yoshida  Theorem).  If 

(a)  A  is  a  closed  linear  operator  on  V  such  that  D(A)  is  dense  In  V, 

(b)  (XI— A)“*  exists  for  some  cjeS?  and  for  every  XeS?  such  that  X  >  w, 

(c)  (|(Xl-ArMl< X  >  w, 

A— “CJ 

then  A  generates  a  strongly  continuous  semigroup  T(t)  with  the  norm  IIT(t)ll  <Me*^ 
The  proof  of  this  theorem  is  given  by  [68 j. 

Definition  3.3.  Considering  the  strongly  continuous  semigroup  T(t)  with 

IIT(t)ll  <Me*‘  M>0, 

if  c*;  =  0,  then  llT(t)  II  <  M  and  T(t)  is  called  an  equi-bounded  semi-group  with  l  >  0. 


Definition  3.4.  For  equi-bounded  semigroup  when  M  =  1, 


-  36- 


IIT{t)ll<  1,  t  >0, 

therefore  T(t)  is  called  contraction  semigroup. 

Contraction  semigroups  are  very  important  in  the  stability  theory  of  semigroups. 
The  contraction  semigroups  are  closely  related  to  the  dissipative  property  of  the 
infinitesimal  generator  A  of  the  semigroup  T.  To  study  the  dissipative  property  of  an 
operator  A,  one  needs  to  apply  inner  products.  Hence,  the  Hilbert  space  would  be  a 
natural  space  for  the  study  of  the  dissipative  property  of  operators. 


Definition  3.5.  An  operator  A  defined  on  D(A)  C  H  (Hilbert  Space)  is  dissipative  if 
Re<Av,v>  <  0  for  every  veD{A). 


Theorem  3.4;  (Phillips  and  Lumer  Theorem).  Letting  A  be  a  linear  operator  with 
domain  D(A)  and  range  of  A,  both  of  which  are  in  the  Hilbert  Space  H  where  D(A)  =  H 
(i.e.,  the  domain  of  A  denses  in  H),  then  A  generates  a  contraction  semigroup  on  H  if 
and  only  if  A  is  dissipative  with  respect  to  the  inner  product  defined  on  H  and  there 
exists  a  X  >  0  such  that  R(XI  —  A)  =  H. 

The  proof  of  this  theorem  is  given  in  [68].  It  should  be  mentioned  that  if  the 
hypothesis  of  the  above  theorem  is  true,  then  one  can  derive 


(XI-A) 


<  X 


-1 


for  all  X  >  0.  Hence,  the  Hille-Yoshida  theorem  could  be  used, 
one  can  find 


By  setting  M  =  1,  =  0 


IIT(t)  II  <  Me^‘  =  1 


which  indicates  that  the  operator  A  generates  a  contraction  semigroup. 


Corollary  3.1. 


If  A  is  a  closed  linear  operator  with  dense  domain  in  H,  then  A  generates 
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a  contraction  semigroup  if  and  only  if  A  and  A*  are  dissipative.  Note  that  A*  is  the 
adjoint  operator  of  A,  given  by: 

<x,Ay>  =  <A*x,y> 

Thus  far,  the  theory  has  covered  the  properties  of  an  operator  A  with  regard  to 
semigroups.  In  general,  the  operator  A  of  an  evolution  equation  can  represent  a 
hyperbolic  or  parabolic  problem.  The  specific  problem  which  will  be  considered  in 
Chapter  5  is  a  hyperbolic  system  of  the  class 

V  =  Lv  ve(L2(n),E“) 

The  states,  v,  form  an  Euchldean  vector  space.  The  vector  valued  function  v,  where  each 
element  of  the  vector  forms  a  function  which  belongs  to  an  infinite  dimensional  Hilbert 
space,  is  defined  on  one  dimensional  spatial  region  (domain)  H  =  x£[0,^j.  The  operator  L 
is  similar  to  the  operator  A  previously  discussed  and  is  unbounded  with  the  specific 
form: 

Lv  =  A  +  Bv 
ox 

It  is  possible  to  establish  the  solution  of  this  system  of  evolution  equations  in  terms 
of  semigroups.  The  solutions  of  this  general  system,  depending  on  the  conditions 
specified  by  dimension  of  E“,  can  result  in  either  groups  or  semigroups  |69]. 

Dynamical  Systems 

The  physical  systems  are  modeled  by  the  evolution  equation  of  integro-differential 
structure.  It  is  often  difficult  to  investigate  properties  of  the  motion  of  a  physical  system 
based  on  its  solution  derived  from  a  specific  condition.  To  avoid  this  obstacle  in  the 
evaluation  of  selected  properties  of  the  system,  such  as  boundedness,  stability  and 
asymptotic  behavior,  an  abstraction  of  the  physical  systems  is  considered.  This 


-  38- 


abstraction  of  mathematical  models  for  the  evolution  of  physical  systems  was  first 
proposed  by  Zubov  [10]  as  dynamical  systems.  The  significance  of  such  systems  is  that 
Lyapunov  stability  theorems  can  be  applied  to  them. 

Definition  3.6.  A  dynamical  system  on  a  metric  space  V  is  a  mapping  u:  R"*"  XV  — +  V 
such  that 

(i)  u(',v)  :  R"^  — ►  V  is  continuous  (right  continuous  at  t  =  0) 

(ii)  u(t,  •)  :  V  — ►  V  is  continuous, 

(iii)  u(o,v)  =  V, 

(iv)  u(t-f4,v)  =  u(t,  u(s,v)), 
for  all  t,  seR'*",  veV. 

For  a  dynamical  system  on  a  metric  space  V,  the  mapping  u(-,v)  :  R"'"  — ►  V  is 
called  the  motion  starting  from  veV. 

From  the  properties  of  a  Co-semigroup  in  definition  3.1,  it  is  clear  that  every  Co- 
semigroup  T(t)  determines  a  dynamical  system  and  vice  versa. 

Definition  3.7:  (Linear  dynamical  systems).  A  dynamical  system  T(t)  on  a  Banach  space 
V  is  linear  if  T(t)€i?(V,V)  for  every  teR'*’,  i.e.,  T(t)  is  a  bounded  linear  operator  on  V. 
From  the  study  of  semigroups  it  can  be  shown  that  for  a  linear  dynamical  system  there 
exists  real  numbers  M,  v,  M  >  1  and  weR  such  that 

llT(t)  II  <  Me^‘  V/teR^ 

A  dynamical  system  may  also  be  defined  on  a  finite  dimensional  state  space  of  a 
finite  order  system.  In  this  case,  the  linear  dynamical  system  is  denoted  by  the 
semigroup  exp(At).  Therefore,  the  generation  of  finite  and  infinite  dimensional 
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dynamical  systems  can  be  investigated  by  using  the  theorems  regarding  the  generation  of 
semigroups  and  from  the  equivalency  between  semigroups  and  dynamical  systems. 

Some  definitions  related  to  the  motion  of  a  dynamical  system  are  provided  below 
to  enhance  the  analysis  of  the  general  motion  of  a  system. 

Definition  3.8.  An  orbit  (positive  orbit)  ^''‘(vo)  in  V  is  the  union  of  all  solutions  of  the 
dynamical  system  T(t)  for  all  times  t  >  o,  i.e., 

7^(v)  =  U  [T(t)Vo] 
t>o 

This  definition  represents  a  specific  trajectory  of  the  system  for  all  times  t  >  o. 

Definition  3.9.  If  T(t)  is  a  dynamical  system  on  V,  then  the  set  11 CV  is  called  a  positive 
invariant  under  T(t),  where  for  every  veil  there  exists  a  T(t)ven  for  all  teR"^.  The  set 
which  is  positive  invariant  under  T(t),  is  called  an  invariant  set.  Clearly,  an 
invariant  set  11  exists  if  and  only  if  the  positive  orbit  7(v)  is  bounded.  In  that  case,  n(v) 
covers  7(v);  therefore,  n(y]  can  be  chosen  as  any  bounded  set  that  covers  7(v).  In  a 
study  of  asymptotic  behavior  of  the  motion  of  a  dynamical  system  as  t  — *-oo  it  is  often  of 
interest  to  find  the  smallest  subset  of  n{v)  such  that  it  has  the  same  properties  as  IT. 
The  task  of  finding  this  smallest  subset  could  be  achieved  by  considering  a  sequence  of 
projections  on  11.  This  leads  to  a  limit  set,  namely,  the  smallest  set  of  invariant  sets. 
Therefore,  this  process  would  become  a  basis  for  the  definition  of  a  positive  limit  set 
below. 

Definition  3.10.  If  for  every  veV  there  exists  a  set  n(v)  such  that  as  — »•  cc,  then 
T(tn)v  — *■  n(v),  that  set  is  called  a  positive  limit  set.  Alternatively,  it  can  be  shown  that 
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n(v)  «  n  [c/  U  T{t)v]  =  n  [C<''y(T(t)v)] 

r<R^  t>r  rtR-^ 

where  C€  mean  the  closure  of  the  set.  Fig.  3.1  shows  the  schematic  of  invariant  sets  and 
limit  sets  for  finite  dimensional  systems  in  V  =  R*  and  an  infinite  dimensional  system  in 


V  s  {[o/],E*). 


*•0 


dlv^,  M)  •  4  d(V,  M)  ■  c 

b)  Invariant  aat  M  in  V  •  E^) 


Figure  3.2 


Invariant  sets  of  finite  and  infinite  dimensional  systems. 
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From  the  definition  of  limit  set  it  is  clear  that  n(v)  CC^'y(v).  mentioned 
earlier,  if  7(v)  is  bounded,  then  n(v)  exists  for  finite  and  infinite  dimensional  systems. 
For  finite  dimensional  dynamical  systems,  the  boundedness  of  7(v)  would  be  sufficient  to 
guarantee  the  existence  of  the  limit  set  n(v)  as  a  limit  set  of  the  invariant  set  n(v). 
However,  for  infinite  dimensional  dynamical  systems,  mere  boundedness  of  ^(v)  would 
not  be  sufficient  for  the  existence  of  n(v).  In  this  case,  the  limit  set  n(v)  is  nonempty  if 
'y(v)  is  compact  (precompact  if  V  is  complete).  More  specifically,  this  can  be  stated  in 
the  form  of  the  following  theorem: 

Theorem  3.5.  If  V  is  complete  and  'y(v)  is  precompact,  then  r2{v)  is  nonempty,  compact, 
connected  and  invariant,  and  T(t)v— ►^{v)  as  t— ►©©.  In  addition,  if  T(t)v — ►HCV  as 
t— *■00,  then  n(v)CC^n  [26].  □ 

Although  most  of  the  aforementioned  materials  can  immediately  be  used  to  cover 
the  invariance  principle,  this  discussion  will  be  continued  after  the  presentation  of  the 
Lyapunov  stability  results.  Some  of  the  preliminary  concepts  and  definitions  which  are 
instrumental  in  Lyapunov’s  direct  method  are  provided  in  the  section  below. 

Definition  3.11.  For  a  dynamical  system  {T(t)}t>o  on  a  metric  space  V,  the  state  Vg  is 
called  the  equilibrium  state  if  T{t)vg  =  Vg  for  all  teR'*’.  In  the  framework  of  limit  sets 
n(vg)  =  {Vg}. 

Definition  3.12.  Letting  {T(t)}t>o  be  a  dynamical  system  on  a  metric  space  V,  the 
continuous  functional  if:V — ►R  is  called  a  Lyapunov  functional  for  {T(t)}t>o  on  V  if 
i?  <  0  for  every  veV.  The  function  if: V  —*0t  is  defined  by 
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k  =  llm  mf  i  (if(T(t)v)  -  if(v)| 

From  this  definition,  it  is  clear  that  ^  does  not  need  to  be  continuous  and  in  this  case,  is 
taken  as  the  right  derivative  of  if. 


Theorem  3.6.  If  {T(t)}t>o  is  a  dynamical  system  on  metric  space  V  and  if:V— is 
Lyapunov  functional  for  all  T(t)veV  for  t£(0,oo),  then 


(i)  if(T(t)v)  is  non-increasing  with 


t 

if(T(t)v)  =  if(v)  -I-  /  i?(T(r)v)dr 


O 


VAe[0,oo) 


(3-19) 


(ii)  In  addition,  if  if(v)  <  —(x5£{v)  for  some  a  >  o,  then 

if(T(t)v)<i^'if(v).  (3-20) 

Proof: 

Part  (i)  -  Since  if(T(t)v)  is  absolutely  continuous  from  the  extended  fundamental 
theorem  of  calculus,  then  equation  (3-19)  will  be  derived  with  its  integration  defined  in 
the  sense  of  Lebesgue.  Therefore, 

t 

if(T(t)v)  -  if(v)  =  /  if(T(r)v)  dr  <  0 

o 

^(T(t)v)  <  i?(,)  (3.21) 

Part  (ii)  - 

t  t 

if(v)  -  i?(T(t)v)  =  /  -i’(T(r)v)  dT>f  a  if(T(r)v)  dr 

O  O 


From  equation  (3-21) 
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^(T(t)v)  <  i?(T(r)v)  V/re[o,tl 

if(v)  -  i?(T(t)v)  >  at  if(T(t)v) 

^(v)  >{l+cA)  if(T{t)v) 

Dividing  t  by  n  and  applying  the  same  inequality  for  duration  of  t/n  results  in 

i?(v)  >  (1  +  a  t/n)  if(T(t/n)v)  >  (1+  cd/nf  if(T(2t/n)v) 

or 

^(v)  >  (1+  at/n)“  if(T{t)v)  . 

As  the  number  of  time  divisions  approaches  infinity,  then  the 

lim  (1  +  at/n)'^  =  e“‘ 

n— *00 

Hence 

-S^’(T(t)v)  if(v).  □ 

If  the  Lyapunov  functional  X  is  defined  on  some  subset  GCY  instead  of  V,  the  subset  G 
must  contain  the  positive  orbit  of  the  system  ^(v);  namely,  G  must  be  an  invariant  set. 

Lyapunov’s  Direct  Method 

The  importance  of  the  Lyapunov  stability  method  when  using  Lyapunov 
functionals  is  to  make  statements  about  the  stability  of  a  dynamical  system.  In  the 
dynamical  system  {T(t)}t>o,  one  must  know  the  characteristics  of  the  mapping  T(t)v  or 
the  trajectory  of  the  system.  The  Lyapunov  approach  is  based  on  making  statements 
about  the  system  characteristic  by  showing  non-increasing  property  of  a  Lyapunov 
functional  if(T(t)v).  This  process  does  not  require  explicit  apriori  knowledge  of  the 
T(t)v.  In  this  section,  preliminary  concepts  and  related  definitions  and  theorems  of 
Lyapunov  stability  are  given. 

Definition  3.13.  For  a  metric  space  V  with  the  metric  d(vi,V2),  where  Vi,V2eV,  the 


distance  of  a  point  from  a  set  KCV  is  defined  as: 

d(vi,K)  =  infd{vi,v)  and  veK  . 

v<K 

If  one  knows  the  '^distance"  between  a  point  and  a  set  in  a  metric  space,  then  the 
stability  in  the  sense  of  Lyapunov  can  be  determined. 

Definition  3.14;  (Stability  of  a  motion).  A  particular  motion  of  a  dynamical  system 
{T(t)}t>o  is  stable  if  for  every  e  >  0  there  exists  <5(e)  >  0  such  that  d{vi,V2)  <  S  implies 
d(T(t)vi,  T(t)v2)  <€  for  every  v^,  V2€V  and  teR"*".  The  motion  T(-)vi  :  R"'"  — is 
called  asymptotically  stable  if  it  is  stable  and  there  exists  <$  >  0,  such  that  d(vi,y2)  <  5 
implies  that  d(T(t)vi,  T(t)v2)— ►©  as  t — ^00  for  every  Vj,  V2eV.  The  motion 
T(’)vi  :  R'*’  — ►  V  is  called  exponentially  stable  if  it  is  stable  and  there  exists  <5  >  0, 
>  0,  and  M(^  <  00,  such  that  d(vi,V2)  <  5  implies 
d(T(t)vi,  T(t)v2)  <  Me~*‘d(vi,V2)  for  every  Vj,  V2eV  and  teR"*". 

Definition  3.15;  (Stability  of  a  set).  A  set  G  in  metric  space  V  of  a  dynamical  system 
{T(t)}t>o  is  stable  if  for  every  e  >  o  there  exists  (5{e)  >  0  such  that  d{v,G)  <  S  implies 
d(T(t)v,G)  <e  for  all  teR'*'  and  veV.  The  set  G  is  called  asymptotically  stable  if  it  is 
stable  and  if  there  exists  5>o  such  that  d(v, G)  <  5  implies  that  d(T{t)v, G)  — >•  0  as 
t  — ►00  for  all  veV.  The  set  G  is  exponentially  stable  if  it  is  stable  and  if  there  exists 
5  >  0  such  that  d(T(t)v,G)  <  Me““*d(v,G)  for  some  >  0  and  M(6)  <  00  and  for  all 
veV. 

Clearly,  if  the  equilibrium  state  Ve  exists  for  a  dynamical  system,  then  the  stability 
of  equilibrium  can  be  regarded  as  the  stability  of  the  set  {(T(t)Ve)}  =  {vg}.  The 
following  theorem  gives  the  stability  conditions  for  equilibrium  based  on  Lyapunov 
functional  approach: 
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Theorem  3.7.  The  equilibrium  state  of  a  dynamical  system  {T(t)}t>o  is  stable  on  the  set 
of  all  states  v  with  d(v,Vg)  <  r  and  r  >  0,  if  there  exists  a  Lyapunov  functional  i?  on  the 
set  G 

G  =  {V  :  d(v,Vg)  <  r}  , 

such  that 

J?(v)>ifK)+f(d(v,v.)) 

The  functional  f(-)  is  a  monotone  function,  f(-)  :  [0,r)  — ►  R'*'  with  f{0)  =  0. 

In  addition  if  .^(v)  <  —  gid(v,Vg)]  for  all  veG  and  some  monotone  function 
s(')  •  g(0)  =  0)  then  the  equilibrium  state  is  asymptotically  stable.  The 

proof  of  this  theorem  is  presented  in  Appendix  A.  □ 

An  extension  of  this  theorem  for  dynamical  systems  can  be  given  accordingly  in 
the  space  of  perturbed  states  from  the  equilibrium,  namely,  veV. 

Theorem  3.8.  An  equilibrium  state  of  a  dynamical  system  is  stable  with  respect  to  a 
metric  d  if  there  exists  a  Lyapunov  functional  if  on  the  set  G  =  {veV  |  d(v,0)  <  r}, 
where  v  =  (v— Vg)eV  (space  of  perturbation)  and  if  satisfies  the  following  conditions: 

(i)  if  is  continuous  and  positive  definite  with  respect  to  d(v,0)  on  the  G  C  V. 

(fi)  is  negative  semidefinite.  Proof  is  by  a  restriction  of  theorem  3.7. 

Theorem  3.9.  An  equilibrium  state  of  a  dynamical  system  is  asymptotically  stable  with 
respect  to  a  metric  d  if  there  exists  a  Lyapunov  functional  X  with  the  conditions  given 
in  theorem  (3.8),  as  well  as  if(d(v,0))  — ►O  as  t  — ►  +oo.  Proof  is  provided  by  theorem 
3.8  and  definition  3.15 
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Lyapunov’s  Method  for  Stability  of  Linear  Systems 

Consider  the  abstract  linear  evolution  equation  of  the  previous  section: 

V  =  Av  v(0)  =  Vo  ,  (3-4) 

where  the  operation  A  generates  a  semigroup  T(t)  on  Banach  Space  V.  One  can  define  a 

stronger  criterion  for  the  stability  of  this  system  by  defining  the  notion  of  exponential 

stability. 

Definition  3.16.  The  system  with  the  above  evolution  equation  (3-4)  is  exponentially 
stable  if  for  a  given  cj  there  exists 

llT(t)ll<Me-‘^S  t>0. 

Note  that  exponential  stability  implies  asymptotic  stability  but  the  converse  is  not 
always  true. 

Theorem  3.10.  Letting  operator  A  in  the  evolution  equation  (3-4)  be  generator  of  a 
semigroup  T(t),  then  the  null  solution  of  (3-4)  is  asymptotically  stable  if  there  exists  a 
Lyapunov  functional  L(v)  such  that  L(v)  >  0  and  L(v)  <  — ^  II  v  iP  for  veD(A). 

As  addressed  in  the  Lax-Milgram  theorem  (2.7),  if  the  II  •  II  is  a  Hilbert  norm  with 
D(A)  =  H,  then  any  Lyapunov  functional  with  bilinear  form  L:  (HXH)  — *■  R  can  be 
made  into  an  inner  product  of  the  form  L  :  (v  X  w)  — ►  <v,  Sw>,  where  S  is  ymmetric 
real  bounded  operator. 

Hence,  an  equivalent  form  can  be  defined  llv  llj  =  <v,  v>2  =  <v,  Sv>  and 
Q!<v,v>  <  llv  ll*  <  /?<v,  v>. 

By  this  equivalence  relationship,  the  existence  of  a  Lyapunov  functional  would 
result  in  the  existence  of  an  operator  S  with  the  aforementioned  properties.  The 
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resulting  Lyapunov  functional,  by  hypothesis,  satisfies 


L{v)  =  L(v,v)  =  <v,  Sv>  . 


L(v)  =  lim 


=  lim 

t— »o 


L(T(t)v,  T(t)v)-L(v,v) 


L((T(t)  +  I)v,  (T(t)  -  I)v 


f-l/t 

l/t  =  2L(v,Av)  =  2<v,SAv>  =  2<v,Av>2 


L(v)  =  2<v,Av>2  <  —O'  llv  iP  <  — -  llv  II2 

a 

In  addition, 

X<V,V>2  —  <Av,V>2  =  <(XI  —  A)v,V>2  <  ll(Xl— A)v  (I2  llv  II2 
or 


<V,V>2 


X  - 


<Av,V>2 


<v,v>2 


<  ii(xi-A)ii2  llv  Hi 


<Av,V>2  ^ 

where  the  hypothesis  — - -  < 


<v,v>2  ~  2a 

<Av,v>2 


lIXI-AlIji  < 


ll(XI-A)-‘  II2  < 


X  - 


<V,V>2 


1 


-l 

* 

< 

X  +  — 

2a 

-1 


X  +  ^ 


2a 


Using  the  Hille-Yoshida  theorem,  this  relation  illustrates  that  operator  A  is  the  generator 
of  a  semigroup  T(t)  such  that 


IIT(t)ll<Me  2“ 

which  is  a  condition  for  exponential  stability. 


The  Invariance  Principle  and  Asymptotic  Behavior 


The  Invariance  Principle  provides  information  about  the  asymptotic  behavior  of 
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motions  as  t  — oo,  whether  or  not  an  equilibrium  exists. 

As  shown  in  theorem  (3-5),  if  the  positive  orbit  "/(v)  is  precompact  in  a  complete 
state  space  V,  there' exists  a  limit  set  n(v)  such  that  the  motion  of  the  dynamical  system 
at  infinity  will  approach  it.  The  method  used  to  locate  this  limit  set  in  the  absence  of 
explicit  knowledge  about  the  motion  of  the  dynamical  system  {T(t)}t>o  is  provided  by 
LaSalle’s  Invariance  Principle  [21|. 


Theorem  3.11:  (Invariance  Principle).  If  {T(t)}t>o  is  a  dynamical  system  on  a  metric 
space  V  and  there  exists  a  continuous  Lyapunov  functional  if  :  V  — ►  on  a  set  G  C  V 
such  that  if  <  — W(v)  for  all  veG,  where  W  :  G  — ►  R'*'  is  a  positive  lower  semi- 
continuous  and  7(v)  C  G,  then  n(vj  C  M'*’.  The  set  M"*"  is  the  largest  invariant  subset  of 
the  set  Ml , 


Ml  =  |v€G  I  if(v)  =  Oj 
In  addition,  if  V  is  complete  and  ^(v)  is  precompact  (compact),  then  T(t)v  — ►  M"^  as 
t  — ►  -j-oo.  As  shown  later  in  this  manuscript,  this  theorem  will  be  used  to  derive 
stabilizability  conditions.  The  proof  of  this  theorem  is  given  in  Appendix  B.  □ 


Stability  from  Spectral  Method 

Amother  alternative  theorem  for  stability  analysis  of  dynamical  systems  originates 
from  the  spectral  properties  of  the  evolution  operators  for  these  systems.  One  of  the 
major  insufficiencies  in  the  application  of  the  spectral  method,  as  opposed  to  Lyapunov 
method,  is  the  restriction  of  the  former  to  linear  dynamical  systems.  For  infinite 
dimensional  systems,  spectral  analysis  gives  the  necessary  conditions,  whereas 
Lyapunov-type  analysis  provides  sufficient  conditions  for  the  stability  of  the  system. 
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For  finite  dimensional  linear  dynamical  systems,  the  spectrum  of  the  evolution 
matrix  A  consists  of  a  finite  number  of  eigenvalues  X  (point  spectrum).  The  system  is 
stable  if  and  only  if  no  eigenvalue  of  the  system  has  a  positive  real  part,  i.e., 

(a)  Sup  Re(X)  <  0  \€a{A) 

The  system  is  asymptotically  stable  if  and  only  if 

(b)  Sup  Re(X)  <  0  \ea{A) 

and  it  is  exponentially  stable  if 

(c)  Sup  Re(X)  <  —S  for  some  S  >  0  and  \€o{A)  . 

For  infinite  dimensional  systems  with  the  evolution  equation  (3-4),  the  spectral  analysis 

would  give  the  following  three  types  of  spectrum: 


(i)  point  spectrum  crp(A),  which  is  the  set 


<Tp(A) 


X  :  (XI— A)  ^  does  not  exist 


(ii)  residual  spectrum  <7'r(A),  which  is  the  set 


<7r(A)  = 


X  :  Range  (XI— A)  is  not  dense  in  state  space  V 


(iii)  continuous  spectrum  <7c(A),  which  is  the  set 


<7c(A)  =  ]  X  ;  Range  (XI— A)  is  not  continuous  |  . 

Requirements  (a),  (b),  (c)  provide  necessary  conditions  for  the  types  of  stability  under 
consideration.  To  relate  these  conditions  to  sufficiency  conditions,  one  can  define: 


=  Sup 


Re  X  :  Xea(A),  a{A)  =  (Tp(A)  n  crR(A)  Pi  o-c(^) 


This  number  is  called  the  lower  index  of  stability.  The  upper  index  of  stability  is 
defined  as  [53] 
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w  =  inf  |gj  ;  llT(t)  II  <  Me-‘ 
where  T{t)  is  the  semigroup  generated  by  the  system  of  equation  (3-4).  In  general,  ^  is 
less  than  cj.  However,  the  condition  where  these  two  indices  are  equal  is  called 
"spectrum-determined  growth  assumption."  Under  this  assumption,  conditions  (a),  (b), 
(c)  will  be  necessary  and  sufficient  conditions  for  their  respective  forms  of  stability.  The 
spectrum-determined  growth  assumption  is  true  when  the  semigroup  T(t)  is  compact 
[54].  It  is  clear  that  the  compactness  of  the  semigroup  T(t)  is  directly  equivalent  to  the 
compactness  of  positive  orbit  of  dynamical  systems.  Therefore,  the  compactness 
condition  is  a  common  necessity  in  both  spectral  and  Lyapunov  stability  analysis. 
However,  spectral  analysis  is  not  directly  applicable  to  nonlinear  systems  unless  the 
system  is  linearized  and  the  analysis  is  performed  in  the  local  sense. 
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CHAPTER  4  -  APPLICATION  OF  LYAPUNOV’S  STABILITY  TO 

MAGNETO-PLASMA  DYNAMIC  SYSTEMS 
Motivation  to  Study  MPD  System 

Recent  increases  in  space  missions  and  the  construction  of  the  space  station  has 
attracted  attention  to  new  alternatives  to  chemical  propulsion  systems.  One  such  system 
is  categorized  and  known  as  the  electrical  propulsion  engine.  In  general,  electric  rockets 
should  be  able  to  develop  considerably  higher  specific  impulses  than  chemical  or  nuclear 
ones.  However,  this  gain  in  specific  impulse  requires  massive  energy  conversion 
mechanisms.  To  avoid  this,  the  electrical  rockets  generally  provide  lower  thrust  for 
navigation  in  low  gravitational  fields. 

The  propellant  of  an  electrical  rocket  consists  of  either  charged  particles, 
accelerated  by  electrostatic  forces,  or  an  electrical  conducting  fluid  (plasma)  accelerated 
by  electromagnetic  and/or  pressure  forces. 

In  the  present  work,  the  modeling  and  analysis  of  steady  magneto-gas-dynamic 
flow  accelerators,  among  other  categories  of  electromagnetic  accelerators,  are  considered. 
The  study  of  MPD  stability  has  been  divided  into  two  parts.  In  the  first,  addressed  in 
the  present  chapter,  the  stability  of  wave  motion  due  to  the  combination  of  an  external 
(applied)  magnetic  field  and  induced  electro-magnetic  field  on  the  plasma,  perturbed 
from  rest  condition,  has  been  investigated.  The  second  part,  addressed  in  the  following 
chapter,  includes  a  study  of  the  stability  of  wave  motion  when  the  flow  is  perturbed 
from  a  non-zero  velocity  equilibrium  state. 

MPD  System  Model 
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The  plasma  dynamic  equations  for  a  magneto-plasma  dynamic  system  are  listed  in 
Appendix  C  in  their  general  form.  These  equations  consist  of  Maxwell’s  equation,  Ohm’s 
law,  conservation  of  electric  charge,  equation  of  state  and  a  set  of  mass,  momentum  and 
energy  equations  [70].  It  has  been  assumed  that  the  plasma  is  originally  at  rest  with 
pressure  Pq,  temperature  Tq,  and  density  with  a  uniform  external  magnetic  field  Hq 
present,  but  no  applied  electrical  field. 

If  the  plasma  is  perturbed  by  a  small  disturbance,  then  the  state  of  the  system  is  a 
combination  of  equilibrium  and  perturbed  states,  hence  the  instantaneous  pressure, 
temperature,  density,  electric  and  magnetic  fields  and  current  density  can  be  written  as 

P=P, +P'(x,t) 

T=To  +T'(x,t) 
p  *Po  +p'(x,t) 

E  =TE,(x,t)  +T Ey(x,t)  H-k  E,(x,t) 

H=T[^H,  -l-hx(x,t)j  -l-T[Hy  +hy{x,t)J  +'kh*{x,t) 

J  =  i  Jx(x,t)  -f  j  Jy(x,t)  k  Jj(x,t)  . 

In  the  case  of  a  neutral  plasma,  i.e.,  cs:  0,  the  number  of  ions  and  electrons  per  volume 
of  plasma  are  nearly  equal.  The  application  of  perturbations  results  in  a  set  of  linearized 
governing  dynamic  equations  of  perturbed  states  [70|.  Following  the  modeling  results 

(9hx  (?hx 

presented  in  Appendix  C,  and  considering  the  fact  that  —  =  0,  it  is  then 

ax  r)t 

possible  to  distinguish  between  two  modes  of  wave  propagation,  the  transverse  mode  (z- 
direction)  and  the  longitudinal  mode  (x-y  direction). 

(i)  Transverse  Mode.  The  equations  governing  the  transverse  mode  are: 
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-  AZ,  Z  -  (h., 


w 


0  <x  <^ ,  t  >0 


(4-1) 


A  = 


<92 

YLa. 

Hx  5x  ^  ax^ 


H 

52 


(4-2) 


1 


V 


where  - - -  and  —  \  X  — Hjj.  The  parameter  V,,  13  defined  as  the  x- 

W  Po 

component  of  AJfven  wave  speed  (70). 

(ii)  Longitudinal  Mode.  The  state  equations  for  this  mode  can  be  reduced  to 


dl 

di 


A  = 


dx. 


vx^  A  A_ 

Hx  ax* 


Vy*  a 


Hy  dx 


4  5* 


dx 

R  a 


=  AZ,  Z  =  j^hy,  V,  u,  /,  t" 


«4 


0  <  X  <  /,  >0 


0 


-RTo 


dx 


K  a* 


Cv  ax 


where  p  =  A,  T  =  _,  Vy  = 
Po  -*•0 


V 


(4-3) 


(4-5) 


'hd 


Hy.  The  parameter  Vy  is  defined  as 


the  y-component  of  the  Alfven  wave  speed. 


Lyapunov’s  Functional  and  Stability  Analysis 
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In  this  section,  the  Lyapunov  functional  approach  is  applied  to  each  mode  of  the 
plasma  dynamics.  The  stability  results  are  derived  and  discussed. 

(i)  For  transverse  mode,  equation  (4-1),  the  boundary  conditions  are 

Z(0,t)  =  0,  Z(^,t)  =  0,  Z(x,0)  =  Zo  6  L2([0,^],E2)  (4-6) 

In  order  to  apply  the  Lyapunov  stability  theorem  to  the  equation  of  motion  for  the 
transverse  mode,  one  should  first  investigate  whether  this  system  generates  a  semigroup; 
namely,  whether  the  wave  propagation  in  the  transverse  direction  represents  a 
dynamical  system.  One  way  to  study  this  property  is  to  use  the  Corollary  of  the 
Phillips-Lumer  theorem,  for  which  the  adjoint  operator  must  be  derived.  The  domain  of 
operator  A  in  (4-2),  D(A)  was  defined  as 

D(A)  »  (h*,w) :  (h*,w)  €  n  X  n  Hi 

V 

The  variables  h,  and  w  each  belong  to  a  set  of  twice  differentiable  functions  in  L2(x),  for 
which  the  first  derivative  in  L2(x)  exists  on  the  boundary  dQ.{x=^,  x=€)  and  dVt  is  a 
compact  support  of  functions  (hj,w).  The  set  of  functions  in  this  domain,  D(A),  are 
dense  in  the  state  space  defined  by  the  norms  Hq  X  Hq-  The  assignment  of  this  normed 
space  is  related  to  the  choice  of  proper  space  for  the  generation  of  dynamical  system  and 
stability  results.  Therefore,  the  inner  product  defined  on  V  =  H  s  XH°  is  <-,->v 


<u,Ay>v.-/|u„u,|  vl 

o  '9(-)  ud^{-) 

(9^vi  +  Hx  dv2 

V* 

77^  (9vi  +  //  <9^V2 

=  i^Ui^Vi  lo  —  i^(9uiVi  lo  +  /  HlViC^^UiClx  +  HxUiV2  Io  ~  /  HxV25ui  dx 

o  O 

'*’  H7  ~  ir“  +  Ml25v2  lo  -  (9u2V2  |o  +  /  Mr2a2u2dx 

(4-8 

After  the  substitution  of  boundary  conditions 


Ui(o)  =  U2(o)  =  Vi{o)  =  V2(o)  =  0 
Ui(^’)  =  U2(0  =  Vi(^’)  =  v^if)  =  0 


The  inner  product  becomes: 


<u,Av>v  =  /  [vj,  vj] 


-^3() 


=  <v,  A*u>v  =  <A*u,  v>v  . 


Hence, 


A*  = 


ui 

U2 


dx 


(4-9) 


If  A  is  dissipative,  then.  Re  <v,Av>  <  0.  If,  instead  of  V  =  H°  XH°,  an  equivalent 
inner  product  space  V  is  considered  as 


-  56- 


=<C.,P.>v 

where  P  is  positive  definite  and  self-adjoint. 


If 


P  = 


0 


0 

% 


O'! ,  aj  >  0 


then 


<z.,Z>v'  =  <Z,PZ>v  (4-10) 

<Z,AZ>v'  =  <Z,PAZ>v 


e 

QTi  i/H  ^(O  Hx  <9(-) 

u 

<Z,AZ>v'  =  /  (hj,w]'^ 

0 

a2 

w 

(4-11) 


-1- w 


v2 

Q!2— 5h,  -I- 


r  dx 


e  ( 

=  Q:ii^h,5hj  lo  -  /  «ii^(5hj^dx  -h  |o  -  /  Q;2i/((9w)^dx 

O  0 


Vx 

o;iHxhj(9w  4-  Q;2— w(9hj 

Wx 


dx 


*  X 

If  and  a2  are  chosen  such  that  aiHx  =a2-— — ,  and  a2  =  1,  then 


a,  = 


Therefore, 
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<L,Al>v  =  -  / 


Hx 


dx 


<  y2 

+  /  a(h*w)  dx 


where  hj(0)  =  w(0)  =  0=w(^)  =  hj(^),  resulting  in 


e 

✓  \ 

V 

2 

<Z,AZ>v'  =  -/ 

0 

Hx 

V  } 

(5hj)^  + 

dx  . 

(4-12) 

It  is  clear  that  operator  A  is  dissipative  in  the  equivalent  state  space  V.  Similarly,  for 
the  adjoint  operator  A*  in  the  space  V, 


e 

v* 

K 

<Z,A*Z>v'  =  /[h^.w]"^ 

W 

0 

— q;2Hx5(*) 

' 

2 

Vx 

Hx 

V  ^ 

(5hJ^  -1-  i/((9w)2 

-/  <9(hxw)  dx 

0  “x 


(4-13) 


or 


( 

(  Tr2 

<Z,A*Z>v'  =  -/ 

0 

Vx 

Hx 

\ 

((9hj2  -f  t<(9w)2 

(4-14) 


This  indicates  that  A*  is  a  dissipative  operator  in  the  inner  product  space  V.  However, 


the  dissipativity  of  A  and  A*  are  found  with  respect  to  the  norm  HqXHo  instead  of 
HoXHq.  To  derive  the  dissipativity  of  A  and  A*  with  respect  to  Hq  XHq,  the  relation 
between  these  norms  will  be  sought.  Let  h^  and  w  be  expanded  in  terms  of  their  Hilbert 


-  58  - 


space  coordinates.  This  expansion  is  possible  because  D(A),  as  defined,  is  compactly 
embedded  in  Hq  XHq  (Sobolev  embedding  theorem  2.9). 


h*  =  S  an(t)Sin 


a— 1 


n— 00 


aTTX 


w  =  E  huWSin 


nTTX 


(4-15) 

(4-16) 


71^ 


o  n— 1  m— 1 
( 


f(dh,f  dx  =  /  S  a,(t)  a„,(t)  Cos-^  Cosi^^  dx 


E  n2^a^(t)/Cos^ii^dx 

n-l  ^  o  ^ 

r  r 

/•  _  2  nTTX  ,  f  2  QTTX  . 

j  Cos^  -  dx  =  J  Sin'^  -  dx 

V  -  C 


/  (^*}*  =  E  a2(t)  /  Sin* -5^  dx 

a-1  ^  o  ^ 


n« 

=  /  S  a*(t)  Sin* 

o  n— I  ^ 

^  ^  71^ 

^  /  E  TT 

o  n— I  ^ 


>/4 


nTTX 


dx 


nTTX 


dx 


E  an{t)  Sin 

nTTX 

dx 

E  ain(t)  Sin 

mTTX 

n-l 

V 

m-1 

dx 


f(dh,f  dx  > 


12  f 
TT 


f  hj  dx 


(4-17) 


The  same  would  be  applied  to  w.  Therefore, 
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<L,M>y  =  -/ 


/  \ 

Vx 

Hx 

V  / 

{ahj*  +  ^aw)2 


dx 


<L,K*l>y/  =  <Z,A1>V  <  J 

^  0 


Vx 

Hx 

2 

hj  + 

dx  (4-18) 


Considering 


=  min  >  0 


(4-19) 


then 


<Z,A*Z>v'  =  <Z,AZ>v' 


(4-20) 


where 


<Z,Z>v' 


e 


/ 


h*  -f  w 


2 


dx 


(4-21) 


From  the  dissipativity  of  A  and  A*  in  V  and  the  fact  that  D(A)  is  dense  in  V,  Corollary 
3.1  leads  one  to  the  conclusion  that  the  operator  A  is  an  infinitesimal  generator  of  a  Co- 
semigroup  in  V.  The  normed  spaces  V  and  V"  are  equivalent,  therefore,  it  is  clear  from 
Theorem  3.10  that  A  also  generates  a  Cp-semigroup  in  V. 


Clearly,  the  Lyapunov  functional  is  the  norm  in  equivalent  normed  space 
namely, 


if  =  <Z,Z>v  =  / 


H, 


h*  +  w2 


dx 


4-22) 


Therefore, 
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k  =  2<Z,Z>v'  =  2<AZ,Z>v' 


(4-23) 


Equation  (4-12)  leads  to  the  conclusion  that  i?  <  0,  i.e.,  the  perturbations  in 
transverse  mode  are  stable.  However,  stronger  results  are  conclusive  from  inequality  (4- 
20),  which  indicates  the  exponential  stability  of  the  system  with 

IIT(t)ll^(,,,)  (4-24) 


where  a/  = 


^min  3-nd 


is  given  by  (4-19). 


(ii)  For  longitudinal  mode  of  wave  propagation,  the  generic  form  of  the  evolution 
equation  (4-3)  is  considered  with  Z  =  ^hy,  v,  u,  p*,  T*j  e  L^(0,^),  Z(0,t)  =  Z(^,t)  =  0. 

The  evolution  operator  is  given  by  (4-5).  An  approach  is  taken  similar  to  the  transverse 
mode  to  construct  the  domain  of  operator  A  and  normed  space  V 


D(A)  = 


Z:  Zi  eH*  nHj  ,  i=l,2,3,4,5 


V  =  h8  X  H3  X  Hg  X  Hg  X  Hg  =  (Hg,E5) 


(4-25) 

(4-26) 


such  that 


<‘.  >V  —  <'/>(Hg,  E’ 


The  adjoint  operator  can  be  found  similar  to  the  case  of  the  transverse  mode  but  with  a 


more  laborious  calculation. 
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vx  d 
H,  ax 

,  a 

Hy  ax 

0 

0 

-h4 

a* 

""ax^ 

0 

0 

0 

A*  = 

0 

4  a^ 
s^'ax^ 

R  a 

Cy  dx. 

0 

0 

+rt4 

0 

0 

0 

0 

0 

K  a^ 

PqCv  dx^ 

(4-27) 


To  show  the  dissipativity  of  A  and  A*,  the  equivalent  normed  space  is  defined  as: 


<\  >V  =  <*,P'>v 


where  P  =  diag(ai,  012,  Q3,  as)  and  aj  >  0, 

H2  tj2  tt2 

=  1>  ^2  =  0(3  =  — ,  —  RT, 

Based  on  these  values  for  ccj’s  and  similar  algebraic  techniques  used  in  the  transverse 
mode,  the  following  can  be  derived: 


^  y 


(4-28) 


<Z,AZ>v  =  <Z,A*Z>v'  =  - 


aii^H 


dx 


+  aoiz 


1^ 

dx 


‘I7 


K 


PoC. 


CUllPH  llhy  11^  -f-  Q^uWvi'F  H - a3i/)lulP  -I-  Qs 


K 


PoC, 


IIT 


(4-29) 

(4-30) 


Therefore,  the  operator  A  generates  a  semigroup  and  perturbations  in  longitudinal 
modes  represent  a  dynamical  system  on  V.  If  one  uses  the  Lyapunov  stability  method 
with 


if  =  <Z,Z>v' 
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if  =  2<Z,AZ>v'  <  0 

which  indicates  the  stability  of  system.  However,  due  to  the  zero  coefficient  for  llp'lF  in 
(4-30),  the  results  cannot  be  extended  to  exponential  stability  and  the  semigroup 
generated  by  A  in  this  mode  is  a  contraction  semigroup,  i.e., 

IIT(x)ll^(,,,)  <1 

Application  of  Spectral  Analysis 

By  means  of  the  separation  of  variables,  the  semi-group  property  generated  by  A 

.  .  .  . . .  fThXj  . 


can  be  constructed.  The  solution  for  Z  =  „  will  exist  if 

TwXw 


th 

x; 

/  Mw  - 

Th 

V2 

x'h 

Xh 

T 

Whenever  T  and  X  independent  functions  of  x  and  t,  respectively,  the  terms  — ,  —  and 
X*’ 

should  be  constant.  This  results  in  X  being  represented  by  a  real  periodic  function, 
x"  X*  iX  T 

Hence,  —  =  ^  ~  application  of  boundary  conditions 

yields  X^  =  utt  for  n  =  ±1,  ±2, .  Therefore,  the  characteristic  equation  for  the 

transverse  mode  can  be  reduced  to 


with  the  following  roots. 
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S 


V 


The  above  expression  for  indicates  that  sup  [Re  a  (A)|  <  0,  which  is  the  necessary 
condition  for  the  equilibrium  solution  of  system  equations  to  be  exponentially  stable,  i.e., 
an  equivalence  to  uniform  asymptotic  stability  for  the  system  of  equation  (4-1). 


The  application  of  the  spectral  analysis  to  the  longitudinal  mode  results  in  the 
following  equation: 


K 


Po  3  Po 


—  + 


Po 


^-bSC 

3  To(7-1)  ^ 


To(^-l) 


x2  \l  , 


+ 


X*  X2 


+ 


SKX* 


To('y-i)  Po 


=  0 


where  Po  is  the  pressure  related  to  Po  and  Tq,  and  7  is  the  specific  heat  ratio  Cp/Cy. 


This  characteristic  equation  does  not  have  a  closed  form  solution.  Therefore,  with 
the  exception  of  very  simplified  and  special  cases  in  general,  the  spectrum  approach  is 
very  complex  and  requires  cumbersome,  symbolic  manipulations  to  provide  stability 
results.  Furthermore,  X’s,  in  general,  can  be  both  positive  and  negative  numbers, 
resulting  in  two  sets  of  characteristic  equations  with  positive  and  negative  coefficients, 
whereas  the  outlined  Lyapunov  approach  would  provide  stability  analysis  for  the  system 
without  the  knowledge  of  system  solutions,  eigenvalues  or  the  form  of  wave  numbers. 


Wave  Speed 


From  the  characteristic  equations  derived  for  transverse  and  longitudinal  modes  of 
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motion  it  is  possible  to  comment  on  the  characteristics  of  the  speed  of  wave  propagation 
in  both  directions.  In  the  case  of  transverse  waves,  the  eigenvalues  are  given  by  the 
following  equation. 


S„ 


In  general,  can  be  written  as  a  combination  of  real  and  imaginary  parts. 


Sq  =  Re{Sn)  +  i  Im(Sn)  . 


This  results  in  a  wave  with  frequency  =  Im(Sn),  where 

4  r* 


Obviously,  in  order  to  have  a  wave  (under  damped  conditions),  u;*  must  be  positive.  If 
cJq  <  0,  there  is  no  wave  propagation  in  the  transverse  mode.  Hence,  wave  speed  can  be 
defined  as 

X2  “  ‘  4 

or 


2 


1 


which  has  an  elliptical  shape,  as  shown  in  Figure  4.1.  When  l\[  =  u,  the  maximum 
value  of  the  transverse  wave  speed  is  limited  by  Alfven  speed  in  the  x-direction. 
Moreover,  for  the  transverse  wave  to  exist,  the  maximum  difference  in  ~  ^')  is  caused 
by  minimum  (i.e.,  X^  =  rr).  Therefore,  if 


Figure  4.1  Wave  s 
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|i>H  —  t/|  < 


TT 


1 


the  generation  of  transverse  wave  is  plausible. 

Due  to  the  complexity  of  the  characteristic  equation,  only  the  special  case  where 
u  =  K  =  0  and  ¥=  0  are  considered  for  longitudinal  waves.  The  characteristic 
equation  for  this  case  would  be: 


S'  +  ^S2  [^RTo  +  (V2  +  V2)l  +  ^RToV2  + 

V 


X2 


+S-^  X’RTo 


=  0 


2c2 


In  the  complex  -  —  ■  plane  for  =  0,  the  following  can  be  obtained: 


S'  +  -^S^  hRTo  +  V^l  +  ^RTo  =  0 

2  if  i - 

=  -  [-(^RTo  +  V2)  ±  V(tRTo  +  V-)2  -  4-/RT0  V2 


where  =  Vj  +  Vy.  For  the  case  of  approaching  infinity  S3,  the  third  root  would 
be  zero  and  —  approaches  — O'RTq,  as  shown  in  Figure  4.2.  Hence  for  ¥=  0,  there 
are  two  wave  speeds,  Vf^st  and  VsIq^,  where 


vL  <  j  [{7RT0  +  V2)  +  V{tRTo  +  V2)2  -  4-;RT„  V2 
V,bw  <  y  [(-.RT,  +  V2)  -  V(''RTo  +V2)2“-4-:RTo 
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CHAPTER  5  -  LYAPUNOV’S  STABILITY  AND  CONTROL  OF  MPD 

THRUSTERS 

Modeling  of  MPD  Thrusters 

In  this  section  the  model  for  MPD  thrusters,  used  as  accelerators,  is  derived.  This 
is  a  simplified  model  of  the  MPD  engine  and  other  researchers  are  studying  the 
derivation  of  a  more  complete  model  for  the  system.  The  task  of  complete  MPD 
modeling  is  beyond  the  scope  of  this  research.  The  simplified  model,  derived  in  this 
section,  is  used  to  study  the  feasibility  of  the  proposed  stability  and  control  analysis  on 
the  MPD  thruster.  Figure  5.1  is  a  schematic  of  such  a  system.  The  flow  of  ionized  gas 
enters  the  thruster  and  is  subjected  to  an  electric  field  E  and  a  magnetic  field  B,  which 
are  perpendicular  to  each  other  and  to  the  gas  velocity.  The  electromagnetic 
acceleration  process  is  an  aggregate  of  effects  from  compressible  gas  dynamics,  ionized 
gas  physics,  electromagnetic  field  theory  and  particle  electrodynamics.  The  individual 
analytic  complexity  of  each  of  these  phenomenon  adds  to  the  level  of  difficulty  in  an 
adequate  theoretical  model  for  this  composite  system.  Analytical  progress  normally 
stems  from  simplified  models  which  preserve  the  essential  physical  aspects  of  a  specific 
situation. 

The  description  of  the  motion  of  plasma  in  terms  of  Maxwell-Boltzmann 
distribution  function  is  too  detailed  to  be  useful  for  many  practical  problems  in  the 
electromagnetic  acceleration  process.  In  these  cases,  the  ion  zed  gas  medium  can  be 
considered  as  a  continuum  fluid  whose  macroscopic  physical  properties  may  be  described 
by  the  conservation  laws  and  Maxwell’s  equations.  These  governing  equations,  for  the 
motion  of  plasma,  will  be  gas  dynamic  equations  with  the  interaction  terms  due  to 
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electromagnetic  forces.  With  this  approach,  a  simplified  model  for  the  problem  can  be 
derived,  as  depicted  in  Figure  5.1.  As  shown,  the  plasma  is  flowing  through  the  constant 
area  channel  along  the  x-axis.  The  channel  is  formed  by  two  conducting  walls  connected 
to  the  cathode  and  anode  poles,  respectively.  Between  these  walls  an  electric  field  E(x)  is 
maintained  in  the  y-direction.  Normal  to  this  electric  field  is  a  magnetic  field  B(x), 
applied  in  the  z-direction.  The  model  is  assumed  to  be  one  dimensional,  i.e.  only 
variations  in  the  x-direction  are  considered.  Applying  this  assumption  in  the  general  set 
of  governing  equations  in  Appendix  C,  one  finds 


Conservation  of  mass: 


^  ,  <9(pu)  ^ 

at  ax 


(5-1) 


where  p  and  u  are  the  density  and  the  velocity  of  the  gas,  respectively.  The  momentum 
equation  results  in: 


<9(pu) 

at 


+ 


a(pu^) 

dx 


--f +F.+F. 


-f 


dT 

ax 


where  p  is  the  thermodynamic  pressure,  r  is  the  shear  stress,  Fg  is  the  electromagnetic 
force  (Lorentz  force)  per  unit  volume,  i.e.  Fg  =  J  X  B,  and  Fy  represents  a  combination 
of  coilisional  forces  and  is  assumed  negligible  compared  to  the  other  terms.  It  is  also 
assumed  that  the  shear  stress  t  is  negligible.  Expansion  of  the  momentum  equation  and 
substitution  from  mass  conservation  equation  results  in  the  following: 


P 


au 

at 


au  ,  m 

-|-  pu-^T - - TT - h  JB 

ax.  ax 


(5-2) 


Due  to  the  one-dimensional  assumptions,  the  variations  in  the  y  and  z  direC^'ons  are 
neglected,  and  the  gas  velocity  components  in  y  and  z  directions  are  omitted.  The 
energy  equation  results  in  the  following: 
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where 


|.(pe)  +  =  JE  + 

di  OX  ox 


e  =  specific  total  energy  (internal  energy  +  potential  energy 
+  kinetic  energy) 

h  =  specific  total  enthalpy  (h  =  e  +  P/p) 

Q  =  heat  flux  due  to  convection  and  radiation,  assumed  negligible. 


On  the  right  hand  side  of  the  energy  equation  above  the  term  JE  represents  the 
Joule  heating  due  to  the  application  of  electric  field.  This  term  can  be  considered  as  the 
dominant  form  of  dissipation  of  energy.  The  plasma  can  be  assumed  as  a  perfect  gas, 
with  the  state  equation 

P  =  pRT,  R  =  R\/m 

where  temperature  is  denoted  by  T,  and  R^  and  m  are  the  gas  constant  of  plasma  and 
molecular  weight  of  the  plasma,  respectively.  As  a  result,  the  specific  energy,  e,  and 
specific  enthalpy,  h,  can  be  represented  as 

e  =  CyT  +  u^/2 
h  =  CpT  +  u^/2 

assuming  the  potential  energy  terms  are  negligible;  Cp  and  Cy  are  the  specific  heat 
coefficients 


Cy  =  R/(7-1) 

Cp  =  '7R/('y-l) 


and  7  is  the  specific  heat  ratio  of  the  plasma. 
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Hence,  the  energy  equation  can  be  simplified  as  follows: 


dt 


u 


PCyT  -1-  P  — 


+ 


dx 


puCpT  +  P 


u 


=  JE 


(5-3) 


From  a  practical  point  of  view,  it  is  reasonable  to  assume  that  the  displaced  current, 
deE 


dt 


-,  and  excess  electric  charge,  p,,  are  negligible  terms  and  that  the  energy  in  the 


electric  field  is  much  smaller  than  that  of  the  magnetic  field.  This  results  in  the 
following  simplified  form  of  Maxwell’s  equations  in  Appendix  C: 


—  •  Vx  B  =  J 

/^e 


where  B  = 


Vx  E  = 


V  •  B  =0 


^B 

dt 


(5-4) 


V  •  E  =0 


If  one  specifies  that  the  magnetic  field  generated  by  the  current  flowing  in  the  gas 
is  negligible  with  respect  to  the  applied  field  B(x),  then  the  electromagnetic  field 
equations  (5-4)  can  be  decoupled  from  the  dynamics  of  the  plasma  [7l]. 

By  Ohm’s  law,  one  can  relate  current  density  J  with  the  applied  fields  as: 

j  =  a(E  —  uB),  cr  =  (7(p,T) 

where  c  is  a  transport  coefficient  and  is  called  the  electrical  conductivity. 

After  substitution  for 

dt 

following  set  of  equations  (5-5)  to  (5-7)  as  the  set  of  dynamical  governing  equations  for 


<3ii 

,  in  the  energy  equation  (5-3),  one  can  arrive  at  the 
dt 
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the  MPD  thruster: 


dt 


dx. 


=  0 


^  ^  ^  _5l  ^ 

dt  ^  dx  p  dx  p 

dT  ,  RT  du  dT  ,  J(E-  uB) 

-  + - ^  +u—  +  - 


dt 


dx 


dx 


PCy 


where  J  is  given  by  Ohm’s  law. 


(5-5) 

(5-6) 

(5-7) 


Control  of  MPD  Thruster  at  Equilibrium  State 

In  the  set  of  non-steady  equations  derived  in  the  previous  section,  i.e.  equation  (5- 
5)  to  (5-7),  the  parameters  B  and  E  can  be  regarded  as  input  functions  to  the  system. 
One  of  the  crucial  questions  about  the  behavior  of  this  system  is  how  to  characterize  the 
relationship  between  system  response  and  those  inputs.  The  general  response 
characterization  of  these  systems,  which  change  with  both  time  t  and  spatial  coordinate 
X,  is  a  complex  problem.  However,  one  can  break  the  problem  into  steps  by  trying  to 
make  assessments  about  the  equilibrium  states  of  the  system  of  the  partial  differential 
equations.  Hence,  one  can  pose  the  question  of  what  choice  of  inputs  would  lead  the 
system  to  a  set  of  equilibrium  states,  and  under  what  conditions  such  controls  would  be 
inplausible.  In  order  to  answer  these  questions  about  the  equilibrium  states  of  the 
system,  one  must  derive  the  equilibrium  set  of  equations  from  the  original  partial 
differential  equations.  If  the  states  of  the  system  are  represented  by  a  vector  valued 
function  v,  as  following: 

P(t,x) 
u(t,x) 

T(t,x) 


V  = 
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then  the  equilibrium  vector  would  be  denoted  by  v,  where 


PeW 

V.  =  |Ue(x) 

Te(x) 


Ui 

■JB 

JE 

Similarly,  the  inputs  (controls)  can  be  defined  by  a  vector  valued  function  U  as: 


U  = 


At  equilibrium,  system  states  reach  their  steady  state  values  and  their  time 

dv^ 

derivatives  would  be  zero,  i.e.  =  0.  Equations  (5-5)  to  (5-7)  can  then  be  simplified 

dt 

to  the  following  equations: 


d  ,  \  n  ^  J 

(peUe)  =0  or  -J—  =  d 


Pe  ,  Ue 


dx  dx  Ue  dx 

dUe  .  R  d(PeTe)  Ui 


u, 


e  —i - — 

dx  Pe 


RTe  dUe 


Cy  dx 


+  Ue 


dx  Pe 

dTe  _  Uz  -  UeUi 

dx 


Pe^v 


(5-8) 

(5-9) 

(5-10) 


Substitution  of  equation  (5-8)  into  (5-9)  results  in  the  following  equations: 

RT„ 


u*  — 


RTp 


du 


dx 


dUe  dTe  Ui 

dx  dx  Pe 


e  dTe  Uj  -  UeUi 

-j-Ue- -  -  - 


dx 


PeCy 
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dUe  •'/UeUi  -  (n'-l)U2 
p,{Ul  -  ^RTe) 

ill  =  -^eUl  +  (Ue^  -  RTe)U2 

dx  peUeCy(Ue  -  iRT^) 


(5-11) 


(5-12) 


If  the  speed  of  sound  ae  is  defined  by  a-l  =  "^RTe,  then  (5-11)  and  (5-12)  can  be  rewritten 
in  vector  differential  equation  form  such  that 


TUe 

-2 

-1  1 

dve  _  1 

Pe 

1 

Pe 

dx  - 

a* 

nl 

n!- 

RTe 

PeCy 

PeUeCy  ! 

-RTe 

1 

7-1 

1 

u  -  pU, 

nl 

Pe 

TUe 

dVe 

dx 

Pe^v 

Pe 

U 


(5-13) 


(5-14) 


Therefore,  for  any  equilibrium  state  it  is  possible  to  arrive  at  a  control  vector  in 
the  local  sense  with  respect  to  the  coordinate  x.  However,  in  transition  from  subsonic 
flow  to  supersonic  flow,  where  u,  =ae,  i.e.  unity  Mach  number,  the  control  inputs  Ui 
and  U2  become  dependent  on  each  other.  In  this  case,  one  should  choose  Uj  and  U2 
according  to  a  certain  relationship  suca  that  the  transition  from  subsonic  to  supersonic 
velocities  and  vise  versa  would  be  plausible,  namely. 


U2 


(7-1)  ’ 


at  Ug  =  ^RTe 


or 


E  _  7Ue 

B  ~  (7-1) 
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This  effect  is  the  same  as  the  choking  condition  in  gas  dynamics,  which  is  extended  to 
magneto  gas  dynamics  [72]. 


Based  on  the  definition  of  Mach  number  as  M  = 


or 


it  can  be  shown  that 


.  .2  ,  dM*  duf  dTe 

M'*  =  —  and 


^RTe 


u 


dM^ 

m2 


=  2- 


du.  dT- 


u. 


Te 


1  d(M2)  ^  ^ 

m2  dx  Ue  dx 


J_  dT^ 

Te  dx 


(5-15) 

(5-16) 


dUe  dTe 

By  substituting  for  — ; —  and  — —  from  equations  (5-11)  and  (5-12),  one  can  derive  the 

dx  dx 


following  results: 


(1-M2)  dM2  ^  ^M2  +  1 
m2  dx  /^e^e^pTe 


(7-1) 


m2  -t-  1 


Pe 


(5-17) 


For  subsonic  flow  M  <  1,  an  increase  in  M  is  possible  provided  the  following 
inequality  for  control  inputs  is  satisfied: 


U2  ^  [(7-1)M2  +  2|  PeUeCpTe  ^  [(7-1)M2  +  2}  ':Ue 
Ui  [7M2  Ij  Pe  [7M2  -h  1]  (',— 1) 

For  supersonic  flow  M  >  1,  an  increase  in  M  along  the  thruster  occurs  if 

U2  ^  [(7-1)M2  ^  2]  TUe 

Ui  [7M*  +  1]  7-1 


(5-18) 


(5-19) 


At 


sonic  condition  the  ratio  of  control  inputs 


U2 

Ui 


should  comply  with  the 
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aforementioned  quantity, 


^/Ue 


7-1 


Clearly,  for  the  decelerating  flow  (when  M  decreases 


along  the  flow)  the  direction  of  inequalities  (5-18)  and  (5-19)  would  be  reversed. 


One  approach  to  constructing  acceptable  control  inputs  Ui  and  U2,  based  on  the 
exclusion  of  singularity  at  choking,  is  to  consider 


[(-(-1)M=  +  2|  I 

^  hM'  +  1|  (-r-l)  '  (1-M*) 


(5-20) 


where  U4  >  0.  Substitution  of  (5-20)  into  (5-11)  and  (5-12)  will  result  in  the  following 
equations: 


dUe  _  UeUi  ^  (7-1)  U4 

dx  pe(7M2-{-l)  TPe 

dT,  _  -2Ui  (M^-1/7)U4 

dx  PeR{7M*-|-l)  PeCyUe 


(5-21) 

(5-22) 


Since  the  Joule  heating  (U2)  is  a  positive  function,  for  supersonic  flow,  from  condition 
(5-18)  one  would  have: 


— 

7UeUl 

,  , 


< 


[(7-1)M^  +  2] 
(M^-l)(7M2-i-l) 


(5-23) 


For  subsonic  flow,  i.e.,  M*  <  1,  U2  in  equation  (5-20)  is  positive  for  any  U4  >  0.  Also, 
the  current  density,  J,  should  be  positive  in  the  chosen  direction  to  accelerate  the  flow; 
hence,  one  can  find  the  following  conditions  from  Ohm’s  law: 


J  =  o{E  —  uB)  >  0 

E  >  uB  (3-24) 


or 
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U2  >  UjUj 


From  equation  (5-20),  U2  can  be  substituted  into  the  inequality  (5-24).  Hence, 


for  M  <  1 


for  M  >  1 


U,(1-M2)  {(^_i)m2-^2)  (7-1)^. 

/  '  «  ~o  ^ 


7u«Ui 

H 

1 

1 

\ 

c 

7-1 

\ 

1 

U4 

^  [(7-1)M^  -I-  2]  (7-1) 

7tleUi 

'  (7MHi)(M2  -  1)  7(M2  -  1) 

7-1 

V  / 

(5-25) 


(5-26) 


It  is  clear  that  in  subsonic  flow,  the  inequality  (5-25)  is  always  satisfied.  Therefore,  there 
is  no  constraint  on  in  this  regime.  However,  for  supersonic  flow,  the  inequality  (5-26) 
represents  a  more  restrictive  constraint  on  U4  than  the  inequality  (5-23).  Therefore,  a 
proper  choice  of  U4  can  be  selected  to  satisfy  the  inequality  (5-26).  The  steady  state 
response  can  be  found  from  equations  (5-21)  and  (5-22)  for  some  arbitrary  choice  of 
Ui  >0  and  by  the  selection  of  U4  according  to  the  aforementioned  process.  It  would  be 
an  interesting  proposition  to  apply  the  theory  of  optimal  control  to  find  the  optimal 
control  inputs  Ui  and  U4  among  an  arbitrary  class  of  functions.  These  optimal  values 
for  Uj  and  U4  can  be  obtained  in  terms  of  the  minimization  of  a  cost  function.  For 
example,  the  cost  function  can  be  selected  from  the  group  of  "fuel  optimal"  problems, 
i.e.. 


x=/r 


IP=  / 


u;  +(u.c„,u,)' 


x=0‘- 


dx 


keeping  in  mind  that  the  optimal  control  problem  is  involved  with  control  constraints  of 
the  form  of  the  inequality  (5-26),  while  U4,  Ui  are  positive  quantities. 


Perturbation  of  Nonlinear  Unsteady  Equations 
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Based  on  the  equilibrium  state  represented  by  Ve(x),  it  would  be  interesting  to 
evaluate  the  system  behavior  in  the  neighborhood  of  its  equilibrium  states.  If  the 
perturbation  of  states  with  respect  to  equilibrium  is  denoted  by  v,  then 


Pe(x) 

V  =  V  +  Ve  = 

u(t,x) 

+ 

Ue(x) 

t{t,x) 

Te(x) 

where 

U  <$C  Ue 

t  «Te 

By  substitution  of  v  into  the  dynamic  equations  (5-5)  to  (5-7),  one  would  obtain  the 
following: 


+  (pue  +  uPe  +  PeUe  +  pu)  =  0 


,  f  ,  -X  5(Ue+u)  R  <9(0eTe  +  pTg  +  PgT  +  p  T) 
—  +  (Ue  +  u)  - — -  -b - - - 

OT  cbc  M  n  ^ 


u, 


Pe  +  P 

df  ,  R(T,  +t)  a(Ue-bu)  ^  ^ 

+ - - h  (Ue  +  u) 


Pe  +  P 

a(Te  +  f )  U2  -  (Ue  +  u)Ui 


dt 


dx 


dx 


(Pe  + 


In  order  to  simplify  the  notation,  perturbation  states  are  denoted  without  the  "  sign 
and  should  not  be  mistaken  for  their  state  functions.  Moreover,  the  nonlinear  terms  are 
assumed  negligible  in  a  sufficiently  small  neighborhood  of  the  equilibrium  states.  .-\_s  an 


example,  one  observes 
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.(9u  dn  dn 


Similar  statements  would  be  true  about  other  nonlinear  terms.  Integration  of  these 
results  with  the  steady  state  equations  (5-8)  and  (5-10)  yields 


dp  ,  dp  .  due 

IT 


+  Pe  -T - 


du  Pe  dUe 


*  cbc  Up  dx 


u  =  0 


(5-27) 


au  ,  RTe  dp  dUe  R  dTg 

“al ^  "j —  P  d"  P 

at  Pe  dx  Pe  dx  pe  dx 


du  dUe  ar  R  dUe 

+  Ue^  +  u+R^  -  iL-_T=0  5-28 

dx  dx  dx  Ue  dx 


ax  ,  ,  ,  Te  dUe  Ue  dT*  Uj 

at  Pe  dx  Pe  dx  PeCv 


RTe  du  ^  dUe  ^  ^  ar  ^ 

+  -  -aT  +  -7—  T+Ue^=0 

Cy  dx  dx  dx 

V  /  V  / 

This  set  of  perturbed  equations  are  linear  and  their  characteristics  depend  on  both 


(5-29) 


equilibrium  states  and  their  derivatives. 


Lyapunov’s  Stability  for  DPS  Applied  to  MPD  Model 

In  this  section,  the  Lyapunov  stability  theorem  is  modified  and  applied  to  the 
system  of  partial  differential  equations  of  (5-27)  to  (5-29).  These  equations  represent  the 
dynamics  of  the  perturbation  state  v  about  an  equilibrium  state  vector  Vg.  Hence,  by 
Lyapunov  method’s  one  can  determine  whether  any  deviation  from  the  equilibrium  state 
is  stable  or  grows  unboundedly  with  time.  To  be  consistent  with  the  notation  in 
Chapter  3,  one  could  rearrange  equations  (5-27)  through  (5-29)  to  the  canonical  form  of 

the  evolution  equation.  Here  the  notation  ( )  is  used  to  represent  —  -  ^ : 


0 
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0 


3v 

at 


Ue 

RT. 


R 


RT, 

0  -  u. 


dy 

dx 


1.  .  R 

—  u*  +  —  T, 
Pt  Pe 


U- 


(7-1)  -  u.  4-  -  T.  — — 

Pt  Pe  PjCy 


+  T, 


(5-30) 


It  is  clear  that  v  =  0,  i.e.,  the  null  solution,  represents  the  condition  where  there  is  no 
variation  from  the  equilibrium  state. 


The  characteristic  directions  (eigenvalues)  of  this  system  of  equations  are: 

^1  =  Ue  .  Xj  =  Ug  -i-  V'^RTe  ,  X3  =  Ue  -  V^R^  > 

Therefore,  (5-30)  represents  a  system  of  hyperbolic  partial  differential  equations.  For  a 
general  hyperbolic  system  of  the  form 


dv 

di 


Adv 

dx 


-f  B  V  L  V 


(5-31) 


where 


dv 

Lv=A  —  +Bv 
dx 


(5-32) 


it  is  possible  to  show  that  the  system  of  equation  (5-30)  can  be  made  into  the  form  of 
(5-31)  with  a  symmetric  A  matrix. 


For  simplicity,  and  are  denoted  by  subscripts  (•)(;  and  (■)x.  As  shown 

dt  dx 

below,  by  dividing  equation  (5-25)  by  (pg),  (5-28)  by  (\/RTe ),  and  (5-29)  by 


(X/ 7—1  Tg),  the  following  symmetric  A  matrix  can  be  obtained: 
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Pt 

p. 

VRTe 
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YrtT- 

Ue 

V('»-1)RT. 

Tt 
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0 
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Ue 
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u 
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T 
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It  can  be  shown  that 


Pe 

U, 

VRTe 

Tx 

Te  V'y-l 

Te  V7-I 


(P/Pe)x 

PPe 

2 

Pi 

\ 

U 

1  1  U  T'e 

Vrt, 

> 

,  2 

+ 


a/  "/—I  T 


(5-33) 


(5-34) 


The  vector  function  v  is  then  redefined  by  following  new  vector 


V 


u 

Vr^ 

T 

TeV-,-1 


(5-35) 


Substitution  of  (5-34)  into  (5-33)  results  in 


Vt-f-Avx+Bv=0 


(5-36) 
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where 


A  = 


u* 


VRT;  Ue  V(^-l)RTe 

0  V(7-l)RTe  Ue 


=  A' 


(5-37) 


and 
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u,  T. 


(7  l)T.u.^u,T',  V(7-1)RT, 

T,V7-1 
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U,  ,  )  1 

.  u.T', 

P,Cv  j 

T.‘ 

T, 

(5-38) 


To  construct  a  Lyapunov  functional  for  this  system,  it  is  sufficient  to  select  the 
functional  as  an  equivalent  norm  of  a  Hilbert  space  of  states  of  equations  (5-36),  i.e. 
from  the  discussion  in  Chapter  3  this  functional  would  be  a  bilinear  form  as  follows 


€ 

^  =  <v,  v>2  =  <v,  Sv>  =  J  v^  S(x)  V  dx 

0 


(5-39) 


where  S  is  a  symmetric  positive  definite  and  bounded  linear  operator.  In  order  to  satisfy 
the  sufficiency  condition  of  Zubov’s  theorem  for  the  stability  of  system  of  (5-31),  one  has 

d^ 

to  show  that  there  exists  a  S(x)  for  which  ~ —  <  0. 

dt 

Equation  (5-39)  yields: 


ASP  ^  T  ^ 

— —  =  /  V  S(x)  V  dx  -f  /  v^  S(x)  V  dx 

Qt  n  n 


Conjugate  operators  are  defined  based  on  the  bilinear  form  operation  as  follows 


I 
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<z,Sy>  =  <S  z,y> 

where  S*  is  (iefin«d  as  the  conjugate  of  S.  Since  S  is  a  real  operator,  S  =  S^,  and 


then 


<  e 

f  y'^  S(x)  V  <ix=  <v,Sv>  =  <S^v,  v>  =  /  (S^v)^  v  dx 
0  0 


=  J  S(x)  V  dx  +  J  (S^(x)  v)"^  V  dx 
^^0  0 

=  <v,  Sv>  +  <S^v,  v> 


If  the  operator  S(x)  is  chosen  to  be  symmetric  S  =  S^,  then 


^  =  2  <Sv,  v>  =  2  f  v’^  S(x)  V  dx 

dt  Q 


(5-40) 


Now,  from  equations  (5-36),  one  can  substitute  for  v  (i.e.,  v.  )  into  (5-40).  This  results  in 


dif 

dt 

di? 

dt 


=  2  <Sv,(-Av,)>  +  2  <Sv,(-Bv)> 

€  e  ' 

=  — 2  J  S(x)  [A  Vj]  dx  —  2  J  S(x)[B  v|  dx 


where 


/  S(x)  A-|^  =  J  ^  [v"^  S(x)  A  v]  dx  -  /  v'" 

avT 


dx  g  dx 


^S(x) 


dx 


A  V  dx  -  J  V  ‘  S(x) 


(:»A 


V  dx 


€ 

-/ 

0 


S(x)  A  V  dx 


(5-41] 


Using  the  same  discussion  about  conjugate  of  the  operator  [SAl,  the  following  can  be 


derived: 
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and 

/FI  €  (  \'^ 

J  S(x)  A  V  dx  =  /  •  (S(x)  A]  V  ■  •  -^  dx  .  (5-43) 

0  [  C7X  J  Q  j  OX 

If  [S(x)  A]  =  [S(x)  A]”^,  equation  (5-42)  becomes  identical  to  (5-43)  and  substitution  into 
(5-41)  results  in 
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CHAPTER  6  -  CONTROLLABILITY  AND  OBSERVABILITY  OF  DPS 

In  the  design  and  analysis  of  a  dynamic  system,  stability  is  often  the  most 
important  objective,  in  terms  of  a  desired  performance.  In  cases  where  the  system  is  not 
stable  or  the  performance  is  not  satisfactory  one  needs  to  know  whether  or  not 
application  of  an  external  input,  namely,  a  control  action,  will  provide  the  desired  result. 
In  case  of  feedback  control,  implementation  of  such  control  action  requires 
observation(s)  of  the  system  behavior.  In  specific  terms,  one  must  be  able  to  determine  a 
control  action  such  that  the  system  can  be  directed  toward  the  desired  performance. 
This  property  depends  on  the  system  characteristics  and  the  way  in  which  it  interacts 
with  the  control  inputs.  On  the  other  hand,  the  acquisition  of  data  may  not  always  lead 
to  the  proper  anticipation  of  the  system’s  behavior.  The  former  property,  which 
determines  whether  a  system  can  be  controlled  to  achieve  a  desired  performance,  is 
called  "controllability".  The  latter  characteristic,  which  determines  whether  the  system 
behavior  can  be  "estimated"  from  observation  of  the  system  output{s),  is  called 
"observability". 

In  this  chapter  these  properties  of  finite  and  infinite  dimensional  systems  and  their 
relations  to  one  another  are  reviewed.  Once  these  properties  of  the  system  are 
investigated,  then  one  can  study  whether  an  unstable  system  can  be  stabilized  by 
addition  of  a  proper  control  action,  namely,  whether  the  system  is  "stabilizable". 

Controllability  of  Dynamical  Systems 

The  abstract  evolution  system  of  Chapter  3  with  the  addition  of  control  action  is 
considered  in  this  chapter.  For  linear  time  invariant  finite  dimensional  systems 
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V  =  Av  +  Bu  (6-1) 

and 

v(o)  =  Vo 

where  vcR®,  ueR“,  Aei?(R‘^,R°)>  and  Be^(R“,R“).  The  solution  of  this  equation  at 
time  t  is 

t 

v(t)  =  e"'-*'Vo  -j-  J  e'^(^“®'Bu(3)  ds  (6-2) 

O 

The  system  of  equation  (6-1)  is  called  controllable  if,  for  each  pair  of  v(t)  and  VoCR"  and 
finite  time  t,  there  exists  a  bounded  control  u  which  directs  toward  v(t).  Without  loss 
of  generality  Vq  can  be  selected  as  zero.  The  set  of  states  v{t)eR“,  which  can  be 
controlled  from  Vq  =  0  by  the  control  input  in  a  finite  time  t,  form  controllability  space 

Theorem  6.1.  The  controllability  space  ^  is  a  linear  subspace  of  R“. 

Proof.  If  Vi  and  V2e'i^CR“,  where  each  state  is  attained  from  the  application  of 
controls  ui  and  U2  at  times  tj  and  t2,  respectively,  then 

Vj  =  J  Bui(s)  ds 

O 

V2  =  J  Bu2(s)  ds- 

O 

Letting  0  <  t^  <  t2,  then  Ui(t)  is  nonzero  for  0  <  t  <  tj  and  zero  elsewhere.  Consider 
an  arbitrary  control  action  u(t)  for  0  <  t  <  t2  such  that 
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u(t)  =  QfUi{t)  +  /3u2(t)  a,^tR  <  oc  . 

Hence,  the  new  control  u(t)  directs  the  system  to  a  new  state  v(t2)e^ 

v(t2)  =  /  e"^**  [<^1(3)  +  '5u2(s)|  ds 

O 

=  j  u,{s)  ds  +  »3v2(t2) 

O 

Letting  s  =  — 12  +  tj  +  s,  then 

s'— t, 

J  Ui(s)  ds  =  J  Ui(3  +  t2— ti)  ds 

where  Ui(s'  +  t2— tj)  in  the  range  of  tj— 12  ^  s'  <  t^  is  the  same  as  Ui(s )  in  the  range  of 

0  <  s'  <  ti .  Therefore, 

tj  S  — tj 

J  gA(t2-S)  ds  =  J  gA.(t,-8)  (jg  =S  V  J  ( t  J  ) 

0  s'-O 

or 

v(t2)  =  QrVx(ti)  + /?V2(t2)  . 

This  states  that  'if  is  a  linear  closed  subset  (subspace)  of  R“.  □ 

A  system  with  state  space  V  =  R“  is  controllable  when  every  point  in  the  space  V 

is  reachable.  This  implies  that  for  every  state  in  space  V  there  exists  a  control  action 
ueU  =  R“.  In  turn,  the  controllability  subspace  must  be  all  of  V  and  the  dimension  of  if 
must  be  equal  to  the  dimension  of  V,  i.e.,  n.  This  idea  can  be  implemented  in  the 
following  theorem  to  derive  controllability  criteria  for  finite  dimensional  state  spaces. 

Theorem  6.2.  The  system  of  equation  (6-1)  is  controllable  if  and  only  if  B*  e^  =  0 
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implies 

V  =  0. 


Proof.  Let  an  input  u{t)  control  the  system  toward  state  Vy  at  time  tj,  where  Vie'€.  If 
state  veV  is  considered  such  that  v  is  perpendicular  to  then 

^1 

<V,Vi>  =  <v,/  e^‘'"“^Bu(s)  ds>  =  0 


=  /  <v,  e^^*'  *^Bu(3)>  ds  =  0 


=  /  <B*  u(3)>  ds  =  0 


Hence,  for  arbitrary  choice  of  ueU,  one  should  have 

B*  =  0  , 


(6-3) 


so  that  the  value  of  inner  product  remains  zero.  For  the  sufficiency  portion  of  the  proof, 
equation  (6-3)  implies  v  =  0.  Therefore,  the  only  vector  perpendicular  to  is  zero,  i.e., 
every  vector  in  V  is  also  contained  in  The  proof  of  necessary  portion  is  obvious, 
because  if  ^  =  V,  v  =  0  and  B*  e^  v  =  0.  □ 

The  well  known  controllability  test  can  be  derived  from  this  theorem.  If  equation 
(6-3)  is  expanded  then 


B*  V  =  B* 


(ti — s)^  , 

V  +  (ti-s)A*v  -I-  - -  A  *2 


2! 


A*  V  -I- 


=  0 


Since  (ti— s)  is  an  arbitrary  time  duration,  then 


l>(ti— s), 


2! 


is  linearly 


independent  vector.  Therefore,  the  coefficient  vector  will  be  zero. 


-  flo¬ 


or, 


B*A* 

B*A*2 


V  =0 


(6-4) 


From  the  Cayley-Hamilton  theorem,  orders  of  A*  from  n  and  higher  are  linear 
combinations  of  orders  (n-1)  and  less.  Therefore,  the  system  of  equation  (6-1)  is 
controllable  if  and  only  if  equation  (6-4)  implies  v  =  0.  This  is  equivalent  to  the 
controllability  test  of 


C* 


B* 

B*A* 

B*A*2 


B*A«-i 


having  rank  n. 


From  the  above  discussion  it  is  clear  that  the  concept  of  controllability  is  basically 
geometric  and  it  is  independent  of  coordinate  system.  To  show  this  let  w  =  Pv  in 
equation  (6-1),  where  P  is  a  non-singular  projection  operator,  then 

P“^w  =  AP“^w  -f-  Bu 
w  =  PAP'^w  +  PBu 


w  =  A  w  -I-  Bu 


C  = 


PB,  PAP-^PB,  (PAp-i)2pB,...,(PAP-^)“-^PB 


C  =P 


B,AB,a2b,...,A“"'B 


-PC 


rank  C  =  rank  PC  =  min(rank  P,  rank  C)  =  min(n,K),  where  rank  C  =  K  and  K  <  n, 
rank  C  =  rank  C.  □ 
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In  the  specific  case  that  P  is  the  similarity  transformation,  PAP  *  =  A 

C  =  I^I,  A,  A^,  •  •  •  ,  PB 

Due  to  the  dlagonality  of  powers  of  A,  matrix  C  has  the  following  form: 

(PB)i  Xi(PB)i  X?(PB)i  Xr'(PB)i 

c  -  i  ;  :  -  : 

(PB)^  XJPB)^  X2(PB)„  Xr'(PB)n 

Clearly,  rank  of  C  is  n  if  and  only  if  there  is  not  any  row  of  PB  which  consists  of  all  zero 
elements 

(PB)i  #  (OlT 

In  infinite  dimensional  systems  the  concept  of  controllability  is  more  complicated 
than  those  of  finite  dimensional  systems.  This  complexity  is  a  natural  outcome  of 
differences  between  finite  and  infinite  dimensional  spaces.  Consider  the  abstract  system 
of (6-1) 

V  =  Av  +  Bu  (6'1) 

v(o)  =  Vo 

where  veV  =  L2([o,^],E“)  and  ueU  =  L2{[o,^],E“)  are  product  infinite  dimensional 
Hilbert  spaces.  Most  of  the  ideas  presented  here  can  be  extended  to  Banach  spaces.  The 
operator.  A  is  assumed  to  be  generator  of  a  semigroup  of  bounded  operators  T(t)  for 
t  >  o  on  the  Hilbert  space  V.  The  operator  B  :  U  — »■  V  is  a  bounded  linear  operator. 
The  solution  of  this  system  as  given  in  Chapter  3  is 
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vi(ti)  =  T(t)  Vo  +  /  T{t^)  Bu(s)  ds  (6-4) 

O 

The  integral  operation  in  (6-4)  is  a  transformation  Cu  of  ueU  into  C  Z.  Equation 
(6-4)  can  be  rewritten  as 

ti 

Vi(ti)  -  T(t)  Vo  =  /  T(t-s)  Bu(s)  ds  =  Cu(u)  (6-5) 

O 

The  left  hand  side  of  (6-5)  can  be  considered  as  a  transformation  of  vg  and  vj  into  Z. 
This  transformation  is  then  S:X  =  (V  X  V)  — ►  Z,  i.e.,  it  maps  the  pair  Vq  and  v(t),  and 

its  image  is  v  =  J  T(t— s)  Bu(s)  ds.  Earlier  in  the  discussion  of  finite  dimensional 
o 

systems,  the  case  of  Vo  =  0  was  considered,  or  in  terms  of  present  notation,  the  space  X 
was  V  X{0}  and  consequently  X  =  V.  Clearly,  there  is  no  loss  of  generality  for  the 
former  representation.  However,  the  present  notation  provides  a  better  insight  in 
distinction  between  Z  and  V,  even  though  they  might  very  well  be  the  same  space. 
Based  on  relationships  between  abstract  spaces  U,X,Z,  the  abstract  linear  control  system 
can  be  established.  Figure  6-1  shows  this  abstract  control  system.  Therefore,  the 
abstract  concept  of  controllability  could  be  addressed  with  the  following  definition: 

Definition  6.1.  The  abstract  linear  control  system  {X,U, Z,Cu,S} 

(a)  is  exactly  controllable  if  ^  =  Range  (€„)  D  Range(S) 

(b)  is  approximately  controllable  if  Range(Cu)  D  Range(S). 

The  definition  6.1  (a)  means  that  for  every  and  Vi(tj)  there  is  a  control  u  such  that 
~  S(vi , Vg)  or  equation  (6-4)  holds.  This  indicates  that  the  control  u  steers  the 
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.1  Abstract  linear  control  system 
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system  from  vq  to  during  [0,ti].  Whenever  Range(S)  is  contained  in  or  is  equal  to  the 
closure  of  Range(Cu),  as  denoted  in  definition  6.1(b),  it  is  possible  that  the  choice  of  vg, 
Vi  would  lead  to  an  image  point  S(vo,vi)  in  Range(S)  which  is  in  Range(Cu)  but  not  in 
Range(C^).  Therefore,  from  relationship  between  Range(Cu)  and  its  closure,  for  a  point 
in  Range(Cu),  there  is  a  point  in  Range(Cu)  which  is  contained  in  the  ball  of  radius  e  for 
every  e  >  0  such  that 

llC,,(u) -S(vo,Vi)  llz  <e  .  (6-6) 

This  shows  that  Range(Cu)  is  dense  with  respect  to  Z.  In  finite  dimensional  system 
Z  =  R*^  both  of  definitions  6-1  (a)  and  (b)  are  equivalent,  since  every  dense  subspace  of  a 
finite  dimensional  vector  space  R“  is  R“.  The  controllability  condition  given  by  theorem 
6.2,  which  indicates  synonymous  approximate  and  exact  controllability  in  finite 
dimensional  systems,  relates  to  only  one  of  the  controllability  concepts  of  definition  6.1 
for  infinite  dimensional  systems.  For  these  systems,  if  there  is  an  element  z  in  Z  such 
that  it  is  perpendicular  to  ^  =  Range(Cuj,  it  can  be  shown  the  z  is  also  perpendicular  to 
Range(Cu),  even  if  Range(C)  is  not  closed.  Consider  an  element  of  Range(Cu), 
defined  as 

Vq  =  J  T(tn-s)  Bu(s)  ds 
0 

Vq  e  Range(Cu)  and  Vn^Range(Cu) 

For  every  e  >  0  there  is  a  V;,  in  Range(Cu)  such  that  Hvjj— v  !l  <  e  where  v  is  the  liniit 
point  of  Vq’s  and  v  ^  Range(Cu).  Then 
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<z,Vtt  >  =  0 

<z,Vn  >  -  <2,v>  +  <z,v>  =  0 

<2,(Vn-v)>  = -<z.v> 

But  v<  =  ("Vn"^)  element  of  Range(Ca)-  Hence,  <z,v>  =  — <Z){vji— v)>  =  0. 

In  the  case  that  e  — 0,  v,  —►0  such  that 

<z,v>  =  lim  <z,v,  >  =  0 
<— 0 

Therefore,  z  is  perpendicular  to  Range(Cu).  Now,  as  in  theorem  6.2, 

t 

<  z,/T(t-s)  Bu(s)  ds  >  =  0 

O 

<B*T*(t-s)  z,  u{s)>  =  0 

where  for  arbitrary  u,  B*T*(t-s)  z  =  0  implies  z  *  0.  In  turn,  this  implies  that  the  only 
element  not  in  the  closure  of  Range(Cu)  is  z  =  0,  namely  Range(Cu)  ^  Range(S). 
Therefore,  the  geometric  concept  used  in  controllability  of  finite  dimensional  systems  is 
only  good  to  derive  the  condition  of  approximate  controllability  of  infinite  dimensional 
systems. 

Observability  of  Dynamical  Systems 

For  the  system  of  equation(s)  6-1,  the  question  of  observability  is  whether  one  can 
reconstruct  the  behavior  of  the  system  from  the  sensory  measurements,  or  more 
specifically,  whether  it  is  possible  to  reconstruct  the  states  of  the  system  from  the 
measured  set  of  outputs.  To  introduce  the  output  set  for  finite  dimensional  system  of 
(6-1),  let’s  consider  the  output  set  defined  by 
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y  =  Cv  (6-8) 

where  ycY  =  R**  and  C  is  a  bounded  operator  from  R“  to  R^,  i.e.,  Ceif(R“,R*^),  where 
p  <  n.  The  output  y  from  (6-2)  and  (6-8)  is 

t 

y(t)  =  -I-  /Ce^*-*)Bu(s)  ds  (6-9) 

O 

since  y(t)  and  u(t)  are  known,  then  y(t)  defined  by 

t 

y(t)  =  y(t)  —  /Ce'^^^“*^Bu(s)  ds  (6-10) 

O 

is  known  and, 

y(t)  =  Ce^‘)v,,  .  (6-11) 

If  the  state  Vq  of  the  system  can  be  determined  from  y(t),  then  the  states  can  be 

reconstructed  from  y(t),  i.e.,  the  system  (6-1)  is  observable.  The  following  definition 
gives  a  more  specific  notion  of  "observability": 

Definition  6.2.  The  system  (6-1)  augmented  with  (6-8)  is  observable  on  [o,ti]  if  given  u 
and  y  as  an  absolutely  continuous  functions,  it  is  possible  to  determine  Vq  uniquely  from 
(6-10).  Clearly,  for  observable  system  of  (6-1)  and  (6-8)  for  an  arbitrary  chosen  u,  the 
output  y  is  zero  if  and  only  if  v<,  =  0.  In  other  words,  the  system  is  observable  if, 

Ce^‘vo  =0  for  o  <  t  <  ti  implies  v^  =  0  .  (6-12) 

Similar  to  the  treatment  of  controllability,  it  is  clear  that  condition  (6-12)  is  the  same  as 
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rank 


C 

CA 


=  Q 


(6-13) 


From  the  controllability  condition  in  theorem  6.2  and  the  observability  condition 
(5-12),  a  duality  relation  between  controllable  and  observable  systems  can  be  concluded. 
This  means  that  if  a  system  is  controllable,  there  exists  a  dual  system  for  it  which  is 
observable.  If  system  (6-1)  is  considered 

V  =  Av  -f  Bu,  v(o)  =  Vq 


The  dual  version  of  this  system  is  its  adjoint.  The  inner  product  <w,v>  between  dual 
states  represents  an  "energy  type"  scalar  quantity.  Therefore,  the  total  energy  of  the 

system  is  / 

0 

performance  index  J  =  <w,v>  I**,  the  following  Hamiltonian  H  would  be  derived: 


dt 


<w,v>dt  =  <w,v>  lo  .  From  calculus  of  variations,  minimizing  the 


H  *  2  Wi  Vi  =  j)  Wi(Av  -h  Bu)i 

i-l  i-I 


Wi  =  =  -^(A*w)i  =>  W  = 


-.\*w 


From  the  transversality  of  the  initial  and  terminal  states  and  since  v  is  not  known 
at  terminal  time  ti,  then  w(ti)  must  be  known,  i.e.,  w(ti)  =wi.  The  only  source  of 
variations  in  energy  function  is  the  control  u  in  original  system  and  the  output  y  in  the 
dual  system.  Therefore, 
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-^<w,v>  =  <y,u> 

<w,v>  +  <w,v>  =  <w,v>  +  <A*w,v)  +  <B*w,u>  =  <y,u> 

Hence,  y  =  B^-w  and  the  dual  of  system  (6-1)  is 

w  =  — A*w  ,  w(ti)  =  y/i  ,  weR“ 
y  =  B*w  yeR“ 

From  the  duality  between  the  observation  and  control,  the  abstract  observed 
system  as  the  dual  to  the  abstract  control  system  of  previous  section  can  be  constructed. 
The  operators  S  and  C^  are  defined  as 

S  :  X=VXV  —  Z 
Cu  :  U  Z 

where  S  and  C^  were  bounded  with  dense  domain.  Therefore,  the  adjoint  operator  of 
those  can  be  defined  as  S*  :  Z  — ♦  X  and  C*^  :  D(C*^)  C  Z  — ►  Y,  with 

<Sx,z>z  =  <x,S*z>x  for  xeX  and  zeZ  (6-15) 

<Cuy,z>7  =  <y,C*ttZ>Y  for  yeD(Cu)  and  zeD(C*u)  (6-16) 

Figure  6.2  shows  this  abstract  linear  observed  system.  This  abstract  system  will  be 
denoted  by  triple  spaces  {X,Y,Z,C*u,S*}. 

The  observability  for  abstract  linear  system  is  more  complex  than  the  specific  case 
of  finite  dimensional  systems.  The  complexity  arises  for  infinite  dimensional  systems  due 
to  the  fact  that  dense  subspaces  of  infinite  dimensional  space  are  not  necessarily  a  whole 
space.  Therefore,  two  concepts  of  observability  would  be  considered. 

Definition  6.3.  The  abstract  linear  observed  system  {X,Y,Z, 0*^,5*}  is 
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(a)  is  distinguishable  if  Ker  C*u  C  Ker  S*, 

(b)  is  observable  if,  in  addition  to  (a)  there  is  a  positive  number  K  such  that 

IIC*nzllY  >  KIIS*zllx  for  every  zeD(C*u)  C  Z  (6-17) 

In  definition  6.3(a),  Ker  C*  is  the  set  of  all  z’s  in  D(C*)  such  that  C*z  =  0. 

The  linear  observed  dual  system  of  an  infinite  dimensional  control  system  (6-1)  is 
w  =  -A*w  =  L2([o,^'],E“) 

y  =  B*w  yeY  =  L2([o,^’l,E“)  (6-18) 

w(ti)  =  Wi 

The  solution  for  y,  as  of  equation  (6-9),  in  the  abstract  sense  becomes 

y  =B*  T*(ti-t)wi  (6-19) 

where  T*(— t)  is  the  semi-group  generated  by  —A*.  The  set  of  all  Wj’s  to  be 
reconstructed,  i.e.,  the  space  of  observed  states,  is  given  by  X  and  the  space  of  outputs  y 
is  Y  and  Z  is  the  set  of  states  of  the  system.  It  is  clear  in  this  case  that  X  =  Z  and  S,S* 
are  identity  on  Z.  The  operator  from  equation  (6-19)  is  found  by 

C\  =B*T*{ti-t). 

In  definition  6.3(a),  the  term  Ker  S*  C  Ker  C*^  means  that  the  set  of  wj  for 
y  =  C*^Wi  =  0  is  contained  in  or  equal  to  the  set  of  all  states  w  =  0.  This  merely  says 
that  if  y  =  0,  then  wj  =  0,  leading  to  the  uniqueness  (distinguishability)  of  observations 
of  states  w.  From  these  relations  it  is  clear  that  the  condition  of  distinguishability  for 
abstract  linear  system  is  that  C*u  Wj  =0  implies  Wj  =  0. 

Definition  6.3(b)  is  much  stronger  than  6.3(a)  since  it  implies  that  for  a  bounded 
output  there  exists  a  bounded  "reconstruction"  of  state,  i.e.,  wjeZ,  or, 
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ltC%  Willy  >KIIWillz 

This  implies  that  there  exists  a  bounded  left  inverse  R  for  such  that 


Therefore, 


R  C*u  Wi  =  Wi  . 


(6-20) 


llwj  ilz  ^  Wi  II 

which  leads  to  (6-20)  with  K  =  The  operator  R  is  the  reconstruction  operator,  as 

llRtl 

shown  in  Figure  8.2.  Therefore,  the  state  Wj  can  be  continuously  reconstructed  from  the 
observation  y  =  C*uWi.  Due  to  this  duality  relationship  between  controllable  and 
observable  abstract  systems,  the  following  theorems  can  be  given: 

Theorem  6.3.  The  abstract  linear  control  system  {X,U,Z,CmS}  is  approximately 
controllable  if  and  only  if  the  abstract  linear  observed  system  {X,Y,Z,C*u,S*}  is 
distinguishable.  The  proof  can  be  seen  from  the  definitions  of  approximate 
controllability  and  distinguishability.  Namely,  lor  every  ueU,  it  was  shown  that  every  z  ♦ 
perpendicular  to  CuU  must  be  zero  to  give  Range(Cu)  D  Range(S).  Hence 
<z,CaU>z  =  0  implies  z  =  0,  but  this  inner  product  is  same  as  <C*uZ,u>u  =  0,  where 
for  arbitrary  u,C*uZ  =0.  Therefore  C*uZ  =0  implies  z  =  0.  Distinguishability  of  the 
dual  system  of  {X,U,Z,Ctt,S}  where  the  observation  operator  is  C*y  is  the  same  as 
C*uZ  =0  implying  z  =  0,  i.e.,  the  controllability  and  observability  is  the  same  for  dual 
systems.  Figure  6.3  shows  this  concept  in  a  diagram  form. 

Theorem  6.4.  The  abstract  linear  control  system  {X,U,Z,Cu,S}  is  exactly  controllable  if 
and  only  if  the  abstract  linear  observed  system  {X,Y,Z,C*^,S*}  is  observable.  The 
proof  of  this  theorem  is  given  in  [40].  □ 
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CHAPTER  7  -  LYAPUNOV-BASED  STABILIZATION  OF  DPS 

One  of  the  interesting  topics  in  control  theory  is  the  relationship  between 
controllability  and  stabilizability  of  dynamical  systems.  When  some  or  all  of  the  states 
of  a  homogeneous  system  is  not  asymptotically  stable,  it  is  desirable  to  control  the 
system  trajectory  such  that  it  becomes  asymptotically  stable,  i.e.,  stabilize  the  system  by 
application  of  a  proper  control  action.  In  this  chapter  stabilizability  of  DPS  is 
investigated.  However,  to  simplify  this  analysis,  the  stabilizability  of  finite  dimensional 
systems  are  first  reformulated  for  ease  of  extension  to  distributed  systems. 

Stabilization  of  Finite  Dimensional  Systems 

Considering  the  finite  dimensional  system  of  Chapter  6  when  the  system  with  null 
control  is  not  asymptotically  stable  and  a  control  action  u  will  be  applied  to  the  system, 
then  the  closed  loop  system  equations  are: 

V  =  Av  +  Bu  veR“,  ueR“  (7-1) 

v(0)  =  Vo 

where  Aei?(R“,R“)  and  B  =  if(R‘“,R“). 

Assuming  u  is  composed  of  the  feedbacks  of  states,  i.e.,  u  =  Kv,  then  this  control 
system  is  stabilizable  if  the  closed  loop  system 

V  =  (A  -t-  BK)v  (7-2) 

v(0)  =  Vo  , 

where  (A  -I-  BK)e  i?(R“,R“),  generates  an  asymptotically  stable  system. 

If  system  (7-1)  is  controllable,  then  it  can  be  shown  that  it  is  stabilizable.  Namely, 
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controllability  is  a  sufficient  condition  for  stabilizability  of  the  system  (7-1).  This  can  be 
shown  by  application  of  controllability  properties.  For  a  controllable  system,  it  was 
shown  that  for  an  arbitrary  u 

B*e'^  1^1  ‘Ij  _  0  implies  2=0 

Moreover,  the  following  inner  product  in  space  Z 

<z,  /  e^^*'“‘>Bu  dt>z  ,  (7-3) 

O 

is  equivalent  to  the  inner  product  in  space  (J 

^1 

/  <B*e^‘^‘‘‘‘^z,u>u  dt  .  (7-4) 

O 

The  corresponding  norm  of  (7-4)  can  be  written  as 

/  <B*e^‘^‘‘~‘^z,  B*e^‘^‘'"‘^z>u  dt  =  /  IB*e^*^‘'"‘^zll?,  dt  (7-5) 

O  O 

For  finite  dimensional  systems,  the  condition  of  approximate  and  exact  controllability 
are  the  same,  i.e., 

lB*e^*(‘'"t)zllu  >5  llzHz  (7-6) 

and  from  equation  (7-5)  it  can  be  shown  that 

f  B*  e^*^~‘^‘'^z,  z>u  dt  >  Jllzllz 

O 

for  some  /?  >  0,  or 
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<  If  (e-^‘B  B*  dt)e^‘‘‘  z,  e^‘*‘  z>u  dt  >  Jllzllz  (7-7) 

O 

An  operator  D  can  be  defined  as 

^1 

D  =  /  B  B*  e-^*‘  dt  ,  (7-8) 

O 

where  D  is  a  self-adjoint  operator.  When  the  system  (7-1)  is  controllable,  inequality  (7- 
6)  and  (7-7)  are  satisfied  and,  consequently,  D  has  a  bounded  inverse.  From  a  Lyapunov 
approach  one  can  show  how  the  operator  D  could  be  used  to  assure  asymptotic  stability 
of  the  closed  loop  system.  In  fact,  at  this  point,  one  can  notice  the  concurrence  of  three 
concepts:  controllability,  stabilizability  and  Lyapunov  stability.  To  show  this,  the 
candidate  gain  operator  for  feedback  is  considered  as  K  =  — B*  Therefore, 

u  =  -B*D'*v  (7-9) 

The  open  loop  system  of  (7-1)  generates  a  solution  of  the  form  v  =  e'^^  v^.  The  closed 
loop  system,  based  on  equation  (7-2)  and  (7-9),  becomes 

v  =  (A-BB*D-‘)v  (7-10) 

v(0)  =  Vo  . 

Consider  a  system  with  the  following  evolution  equation 

y  =  (A  -  B  B*  D-^)‘y  (7-11) 

y(o)  =  Vo 

where  y  is  in  the  space  V  =  R“,  i.e.,  system  (7-10)  and  (7-11)  share  a  common  space. 
The  open  loop  system  (7-1)  generates  the  semigroup  T(t)  and  the  closed  loop  system  (7- 
9)  with  evolution  operator  (A  —  B  B*  D“*)  generates  S(t).  Therefore,  the  system  of  (7- 


-  106  - 


11)  with  evolution  operator  (A  —  B  B*  D  *)*  is  a  generator  of  semigroup  S*(t). 
Lyapunov  function  be  selected  as  the  following  weighted  inner  product  on  VXV 

^{y)  =  <y,Dy>v 
^1 

=  <y,(/  e-^‘  B  B*  e'^'^  dt)  y> 

O 

=  /  <B*e-^*V  ,  B*e-^‘V>u  dt  >  ti  ^Ib^lly  . 

O 

Since  D  is  self-adjoint,  then  the  time  derivative  of  i?  becomes 

^(y)  =  2<y,Dy>v 

=  2<(A-BB’*  D-»)*y,Dy>v 
=  2<A*y,Dy>v  -  2<I>~^B  B*y,Dy>v 

i?(y)  =  2<D  A*y,y>v  -  2lB*yl^- 

where 

D  A*y  =  (/  e-^‘  BB*  e-^*‘  A*  dt)y 

O 

=  -(/  e-At  B  B*^(e-A‘‘)  dt] 

O  dt 

Operators  B  and  B*  are  time  invariant,  hence, 

j 

D  A*y  =  -/  e-^‘  B^(B*  e-^*‘)  dt  , 


Let  the 

(7-12) 


(7-13) 


(7-14) 


and 
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*  A 

<DA*y,y>v  =  -/  B*e'^*‘y>u  dt 

j 

=  — <B*e“^*V.  B*e“^*V>u  fit 

^  o  ^ 

=  -i  lB*e-^*‘'ylb  +  J  IB*ylb  (7-15) 

Therefore,  by  substitution  of  (7-15)  into  (7-14),  it  can  be  seen  that 

i>(y)  =  -iB*e~^‘‘'yllu  -  IB*yllu  <  0  (7-16) 

From  (7-6)  it  is  clear  that 

IB*e"'^‘‘'yllu  >  ^illyllv  (7-17) 

IB*yllN  >(52llyllv  (7-18) 

Therefore  i?(y)  is  not  just  negative  semidefinite,  but  for  finite  dimensional  systems  due 
to  (7-17)  and  (7-18),  it  satisfies  the  more  restrictive  condition  of  being  negative  definite. 

i?(y)<-(^i  +<^2)  Ibrib  (7-19) 

From  this  Lyapunov  functional  and  the  results  of  the  Lyapunov’s  direct  method  applied 
to  linear  systems,  as  shown  in  Chapter  3,  system  (7-11)  can  be  found  exponentially 
stable  with  the 

for  some  a  >  0  (7-20) 

Therefore  system  (7-10)  which  is  the  closed  loop  system  with  the  control  law  of  equation 
(7-9)  is  exponentially  stable,  i.e., 

||g(A-BB*D-')t||  ^  |yA*-D-‘BB*)t|| 


or  from  (7-20) 
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llg(A-BB*D  ‘)t|[  <  e  for  some  a  >  0  .  (7-21) 

Therefore,  controllability  of  a  finite  dimensional  system  leads  to  stabilizability  of  that 
system.  However,  this  relationship  does  not  exist  in  the  case  of  infinite  dimensional 
systems,  as  will  be  discussed  later. 

When  the  condition  for  existence  of  feedback  gain  operator  K  is  given,  it  is  often 
required  that  the  system  characteristics  behave  in  a  desired  manner,  namely,  the 
eigenvalues  of  the  closed  loop  system  being  specified  apriori.  The  theorem  below  proves 
the  existence  of  an  operator  K  such  that  the  closed  loop  system  has  this  placement 
property. 

Theorem  7.1.  Letting  the  system  (7-1)  be  controllable  and  operator  B=b  be  an  nXl 
matrix  (single  input  system),  then  the  closed  loop  system 

V  =  {A+bK)v  (7-22) 

V  =  Vo 

has  the  eigenvalue  placement  property.  This  property  can  be  stated  as:  given  anj  n- 
tuple  of  eigenvalues  (Xi,  X2,  Xq},  there  is  one  and  only  one  row  vect  )r 
K  =  [Ki,  KnI  such  that  the  matrix  A-fbK  has  the  eigenvalues  aj,  X2,  ...,  X^. 

Clearly,  this  theorem  can  be  extended  to  the  case  of  a  general  nXm  B-matrix.  The 
proof  of  theorem  7.1  is  given  in  [73]. 

In  general,  stabilizability  of  system  (7-1)  does  not  guarantee  controllability  of  this 
system.  The  counter  example  would  be  a  homogeneous  stable  system,  which  is  clearly 
uncontrollable.  However,  the  pole  placement  property,  i.e.,  being  able  to  select  the 
eigenvalues  of  the  closed  loop  matrix  A-fBK,  implies  controllability. 
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The  pole  placement  property  of  controllable  systems  can  be  considered  as  an 
algebraic  property  of  matrix  operators.  From  this  property  it  was  shown  that  for  a  set 
of  selected  eigenvalues  of  system  (7-2)  there  exists  at  least  a  matrix  K,  provided  the 
system  is  controllable.  This  means  that  the  characteristic  equation 

det  [XI  -  (A-I-BK)]  =  0  ,  (7-23) 

for  the  given  Xj’s,  i  =  l,2,...,n,  would  result  in  a  set  of  n  characteristic  equations  which 
lead  to  determination  of  K,  provided  system  (7-2)  is  controllable.  Therefore,  the 
controllability  of  the  pair  (A,B)  is  the  necessary  and  sufficient  condition  for  matrix 
(A-l-BK)  to  have  eigenvalues  as  specified  and  this  property  is  independent  of  (A-t-BK) 
being  an  evolution  matrix.  In  the  literature,  the  pole  placement  of  general  matrix 
(A-hBK)  is  not  usually  looked  upon  in  the  framework  of  controllability  of  the  pair  (A,B). 
Despite  the  fact  that  placement  property  of  eigenvalues  of  matrix  (A-I-BK)  and  evolution 
system  (7-2)  is  the  same,  here  this  distinction  is  made  because  of  the  usage  of  this 
property  applied  to  the  matrix  (A-fBK)  in  derivation  of  stabilization  of  a  class  of  DPS. 
These  views  of  stabilization  in  finite  dimensional  systems  set  a  basis  for  their  comparison 
with  infinite  dimensional  systems. 

Stabilization  of  Infinite  Dimensional  Systems 

If  the  abstract  evolution  equation  of  an  infinite  dimensional  system  given  in  the 
form  of  equation  (3-4)  with  evolution  operator  A  generates  a  semigroup  T(t)  which  is  not 
asymptotically  stable,  then  a  control  u  will  be  applied  to  the  system  and  the 
nonhomogeneous  system  can  be  represented  by 


V  =  Av  -I-  Bu 


(7-24) 


-  no¬ 


where 

veV=L2(10,/’],E“) 
ueU  =  L2([0,^’],E“) 
v(0)  =  Vo 

The  state  v  and  control  u  are  vector  functions  of  spatial  domain  xe[0,^].  Hence,  they 
belong  to  n  and  m-tuple  function  spaces  V  and  U,  respectively. 

The  control  system  is  then  stabilizable  if  there  exists  an  operator  K  :  V  — ►  U  such 
that  the  control  u  =  Kv  generates  the  closed  loop  system 

V  =  (A  +  BK)v  (7-25) 

v{0)  =  Vo 

vcV 

and  the  closed  loop  evolution  operator  (A  +  BK)  is  the  generator  of  an  asymptotically 
stable  semigroup  S(t). 

Unlike  the  finite  dimensional  systems,  approximate  controllability  of  system  (7-24) 
does  not  guarantee  stabilizability.  However,  if  the  system  is  exactly  controllable  as 
discussed  in  Chapter  6,  then 

IIC*u  zllu  >  ^llzllz  •  (7-26) 

For  the  linear  system  (7-24),  this  condition  can  be  represented  as 

IB*  T*(ti— t)zliu  >  6  Ilzll2  ,  (7-27) 

for  every  ticR'*'  and  0  <  t  <  tj.  From  this  condition  the  norm 


-  Ill . 


/  zll?j  dt 


will  satisfy  the  following 


(7-28) 


^1 


/ 


l[B*T*(ti-t)  zll?j  dt  llzlll 


(7-29) 


for  some  0  >0. 

The  following  can  be  written  from  the  properties  of  semigroups  T  and  T*. 

T(ti-t)=T(ti)T(-t)  =  T(-t)T(t,) 

T*(ti-t)  =  T*(ti)  T*(-t)  =  T*(-t)  T*(ti) 


Therefore, 


/IB*T*(ti-t)zll?jdt  =  /<B*T*(~t)  T*{ti)z,  B*T*(-t)  T*(ti)z>udt 

O  O 

Similar  to  the  treatment  of  finite  dimensional  systems,  it  has  been  determined  that 

ti  ti 

/IB*T*(ti-t)zll?jdt  =  <(/T(-t)BB*T*(-t)dt)Tiz,  t|z>u  •  (7-30) 

0  O 


If  an  operator  D  is  considered  as 


D(')  -  (/T(-t)BB*T*(-t)dt)  {■) 

o 


then  D(*)  is  a  self-adjoint  operator  and  from  (7-29)  and  (7-30) 

<D  T*(t,)z,  T*(ti)z>u  >  3\\z\\l  .  (7-31) 

From  (7-31)  it  follows  that  D(‘)  will  have  a  bounded  inverse. 
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If  the  control  action  is  considered  as  a  feedback  control  with  the  gain  operator  K 

K(-)  =  (7-32) 

then 

u  =Kv  =  -B*D-‘(v)  (7-33) 

Since  both  operators  D(’)  and  D“^(')  operate  on  a  distributed  state  function  and 
generate  a  distributed  function,  then  u  is  a  distributed  feedback  control.  The  closed 
loop  operator  of  system  (7-25)  generates  a  semigroup  S{t)  and,  similar  to  finite 
dimensional  case,  a  Lyapunov  functional  X  can  be  defined  as 

i?(y)  *  <y.Dy>  (7-34) 

where  y  is  the  state  of  the  following  system 

y  *=(A-BB*D“i)‘y  (7-35) 

y(0)  =  Vo 

yeV  =  L2([0, <’],£“) 


The  closed  loop  operators  of  systems  (7-25)  and  (7-35)  are  adjoint  of  each  other 
and  (7-35)  generates  the  adjoint  semigroup  S*(t).  The  time  derivative  of  if  becomes, 

i’(y)  =  2<DA*y,y>v  -  2<B*y,  B*y>u 

where 


DA*y  =  /  T(-t)BB*T*(-t)A*y 


Since  the  open  loop  evolution  operator  of  system  (7-35)  is  A*,  which  generates  the 
semigroup  T*(t),  where 
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-iT*(-t)  =  -T*(-t)A*  , 

at 

then 

j 

DA*y  =  -/  T(-t)B^(B*T*(-t))  dt 

Finally,  similar  to  the  derivations  of  finite  dimensional  systems,  the  following  can  be 
obtained 

if(y)  =  _|[B*T*(-ti)yll?,  -  llB*yll?j  <  0  (7-36) 

Due  to  the  exact  controllability  and  (7-27),  the  right  hand  side  of  (7-36)  is  bounded  from 
below  in  the  state  space  V,  i.e.,  i?  is  negative  definite  with  respect  to  llyll^.  Hence,  the 
system  (7-35)  is  asymptotically  stable.  It  can  be  seen  that 

IIS*(t)ll  <  Me”*^*  for  some  M  >  1  and  a  >  0 

Since  IIS(t)ll  =  IIS*(t)ll,  then  the  same  statement  applies  to  the  system  (7-25). 

Therefore,  the  exact  controllability  is  a  sufficient  condition  for  stabilizability  of  a 
general  distributed  parameter  system.  However,  the  exact  controllability  condition,  as 
used  in  (7-27),  is  often  very  difficult  to  satisfy  for  general  distributed  systems. 
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class  of  systems,  with  Ai  being  symmetric,  shares  a  common  form  with  the  linearized 
model  of  the  MPD  system  studied  in  this  chapter.  Let  Ai(x)  possess  eigenvalues  of  the 
general  form:  Xj(x)  <  X2W  <  —  <  Xp(x)  <  0  <  Xp+i(x)  <  ...  <  Xij(x).  The 

eigenvalues  show  the  directions  of  the  characteristic  lines  of  partial  differential  equations. 
In  addition,  there  exists  a  continuously  differentiable  matrix  0(x),  based  on  the  system 
eigenvectors  such  that  0“^(x)  A(x)  0(x)  =  A(x).  If  one  considers  a  new  set  of  states  w, 
such  that 


V  =  0(x)w  , 


then  the  substitution  of  v  into  (7-37)  results  in 


+  Ai(x)0(x)-^  + 


A,(x)i^  +B,0(x) 


w  =  0 


■^  +  O  ‘(x)A,(x)0(x)-^  +  0-‘ 


A,(x).^^  +B,0(x) 


w  =  0 


Therefore,  one  can  conclude 


^  +  A(x)f.+/?(x),=0 


(7-38) 


Since  A(x)  is  a  diagonal  matrix  of  the  system  eigenvalues,  and  since  Xj  to  Xp  are 
negative  and  Xp^j  to  X^  are  positive,  then  A  can  be  decomposed  as 


A  = 


A-  O 
0  A+ 


A  =  diag(Xi,  .  .  .  ,Xp) 
A"^  =  diag(Xp+i,  .  .  .  ,XJ 


(7-39) 


The  corresponding  decomposition  can  be  applied  to  the  states 


w  = 


w 


w 


(7-40) 
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Therefore,  (7-38)  can  be  reduced  to  a  set  of  ordinary  differential  equations  of  the 
following  form: 


where  wjj  is  the  kth  element  of  the  vector  valued  function  w.  Moreover,  represents 

the  direction  of  the  kth  characteristic  line,  hence  is  the  directional  derivative  along 

dt 

the  corresponding  characteristic  line.  On  the  right  hand  side  of  (7-41),  0^  is  the  kth  row 
of  the  matrix  valued  function  0{x).  The  set  of  n  ordinary  differential  equations  (7-41) 
are  coupled  by  the  term  Initial  values  for  t  —  0  can  be  written  as 


w(x,0)  =  Wo  €L2([0, /];£“) 


Considering  the  boundary  conditions  at  a  point  on  the  boundary  x  =  0  and  t  =  to , 
as  shown  in  Figure  7.1,  one  finds  the  characteristics  with  negative  and  positive 
eigenvalues  arrive  at  that  point  with  negative  and  positive  slopes,  respectively  [47].  The 
"incoming  information"  consists  of  values  of  wj^  associated  with  characteristic  line 
Ck(O.to)  with  negative  slope  for  k  =  l,...,p.  The  "outgoing  information"  consists  of 
values  of  wj^  associated  with  the  positive  slope  characteristic  Cjj(o,  to),  k  =  p4-l,...,n. 
Hence,  along  the  boundary  x  =  0,  the  values  of  w'*'  should  be  known.  It  is  clear  that 
along  the  boundary  x  =  €,  the  orientation  of  characteristics  will  be  reversed  and  hence 
the  values  of  w"  should  be  known.  When  system  (7-37)  is  unstable  it  can  become 
subjected  to  stabilization  with  addition  of  a  control  action  u,  such  that  the  resulting 
nonhomogeneous  system  will  have  the  generic  form  of  the  abstract  system  (7-24) 


Figure  7.1  Characteristics  configuration  at  x  ==  0  and  t 
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dv 

^=Av+Bu 

dv 

where  Av  »  — Aj  — - Bj  v 

ax. 

v(0)=v„ 

The  input  operator  B  depends  on  the  way  that  the  control  is  interacting  with  the 
system. 


In  general,  the  system  of  equation  (7-37)  is  assumed  to  be  stabilizable  by  some 
form  of  feedback  control  as  u  =  Kv,  where  K  :  v  — *•  u.  Hence,  the  operator  A-l-BK 
would  be  an  infinitesimal  generator  of  an  asymptotically  stable  semigroup.  If  one  defines 
a  Lyapunov  function  as  i?  =  <v,v>,  similar  to  (5-39)  with  S  =  I,  then  the  following 
would  result: 


k  —  — v^AiV  I  +  / 


3Ai 

dx 


-  2Bi  +  2BK 


vdx 


(7-12) 


From  transformation  v  =  0(x)  w  we  have 


v^AiV  I  =w^O^AiOw  I  =w^Aw  | 

O  0  0 


(7-43) 


where 


w 

T 

A- 

0 

w 

w-^ 

0 

A+ 

w"^ 

. 

w^  .\,w  = 


(  e  e 

w^Aw  I  =  w~  A“w~  I  -I- w'''^A''’w'''  I 
0  0  0 


(7-44) 


At  X  =  0,  and  at  x  =  £,  w'^  and  w  ,  respectively,  must  be  known  or  given  from  the 
boundary  conditions,  i.e., 
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(7-45) 

(7-46) 


w^(o,t)  =  DoW~(o,t) 

Hence 

w”^Aw  I  =  w‘^^(^,t)A'^(^)w'^(^,t)  +  w‘*’^(^,t)  D?'A"(^)  w'^(^,t) 

O 

-[w~^(o,t)A~(o)w"(o,t)  4-  w~^(o,t)DjA'^(o)DoW(o,t)] 

To  guarantee  stability  in  the  sense  of  Lyapunov,  the  following  must  be  satisfied: 

-w^''(/,t)  [A+{^')  +  d7A-(0D/)  w^(^’,t)  <  0  (7-47) 

w~^(o,t)  [A“{o)  -h  DJA‘^(o)Do]  w“(o,t)  <  0  (7-48) 

In  the  case  that  w'''(o,t)  =  0  and  w~{^,t)  =  0,  (7-43)  reduces  to 

f 

w'^Aw  I  =w+{/,t)’rA+(0w+(/’,t) 

O 

— (w~(o,t)’^A"(o)w"(o,t)]  >  0  (7-49) 

Hence,  — Ai  v  |  <  0. 

O 

Another  condition  to  be  satisfied  is  that  the  following  matrix  should  be  negative 
definite 

:  negative  definite  (7-50) 

Therefore,  by  using  a  proper  feedback  gain  K,  such  that  conditions  (7-49)  and  (7-50)  are 
satisfied,  the  Lyapunov  functional  leads  to 

if  <  <v,McV>  < -^  Itvllv  (7-51) 

which  results  in  the  asymptotic  stability  of  system  (7-37).  If  Mg  is  negative  definite,  then 
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9Ai 

the  supremum  of  its  eigenvalues  must  be  less  than  zero.  Since  the  operators  — - —  and 

cic 

Bi  are  only  matrix  functions  of  x,  stabilization  of  the  infinite  dimensional  system  ("-37) 
reduces  to  the  pole  placement  property  of  a  finite  dimensional  matrix.  As  mentioned 
earlier  in  the  discussion  of  finite  dimensional  systems,  if  the  pair  (A,B)  is  controllable, 
then  there  exists  a  K  such  that  the  eigenvalues  of  A  -t-  BK  can  be  placed  in  desired 
locations.  Therefore,  there  will  exist  at  least  a  matrix  K  such  that  Mg  is  negative 
definite  if  and  only  if  the  controllability  matrix  C,  below. 


2B  ,  2 

<9Ai 

dx 

-  2Bi 

B,  2 

dAi 

dx 

-2Bi 

2 

B  ,  ....  ,  2 

dAi 

ax 

* 

n— 1 

B 

* 

• 

. 

(7-52) 

has  a  full  rank.  In  the  case  of  the  MPD  thruster,  where  Aj  is  given  by  (5-32)  and  Bj  is 
given  by  (5-38),  its  control  is  a  distributed  (body  force)  type  input.  Addition  of  the 
distributed  control  inputs  can  be  achieved  by  imposing  a  perturbation  to  the  equilibrium 
controls  in  equations  (5-27)  to  (5-29).  Hence,  the  control  vector  in  (7-24)  will  be  defined 
as 


u  = 


Ul 

U2 


(7-53) 


Therefore,  the  input  operator  B,  which  couples  u  to  the  system  of  (7-24),  will  be 


calculated  as 
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B  = 


0 

1 


,Vrt7 


— u. 


0 

0 


PeTeV-/-lC,  PeTeV-^-lC, 


(7-54) 


From  the  values  of  Aj  and  Bi  given  by  equations  (5-37)  and  (5-38),  respectively,  and  B 
given  by  (7-54),  one  can  construct  the  controllability  matrix  for  the  MPD  model.  A 

Te  Ug 

detailed  computation  leads  to  the  condition  that  and  —  should  not  be  zero 

simultaneously  so  that  the  system  of  the  MPD  model  would  be  stabilizable.  This 
indicates  that  although  there  is  not  a  formal  input  term  influencing  the  density  evolution 
equation  (5-27),  i.e.,  the  continuity  equation,  it  is  possible  to  control  and  stabilize  the 
density,  velocity  and  temperature  of  the  plasma  with  the  input  controls  appearing  only 
in  the  momentum  and  energy  equations. 


Stabilization  of  DPS  Represented  by  Contraction  Semigroups 

An  important  class  of  distributed  parameter  systems  which  covers  a  wide  range  of 
hypo-eliptic  and  structural  systems  are  discussed  in  this  section.  Consider  a 
nonhomogeneous  infinite  dimensional  system  represented  by  the  evolution  equation  (7- 
24),  where  the  open  loop  evolution  operator  A  generates  a  contraction  semigroup  T(t), 
defined  in  Chapter  3  as  irr(t)ll  <  1.  From  the  Hille-Yoshida  theorem  it  has  been 
determined  that  the  generator  A  must  be  dissipative  with  respect  to  an  appropriate 
norm  so  that  this  contraction  exists.  For  those  dissipative  generators  with  iT(t)Ii  <  1 
the  evolution  process  is  already  asymptotically  stable  and  stabilization  has  no  meaning. 
However,  for  the  case  where  <v,Av>v  =  0  and  HT(t)ll  =  1,  the  dynamical  system  is  not 
asymptotically  stable  and  stabilization  can  be  addressed.  This  category  of  problems 
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contains  a  majority  of  structural  systems  with  negligible  damping  which  are  often 
subject  to  undamped  oscillations.  Consider  the  system  (7-24) 

V  =  Av  -t-  Bu 
v(0)  =  Vo 
veV 

which  is  subjected  to  the  following  condition:  <Av,v>v  =  0. 


Assuming  a  Lyapunov  functional  if  on  a  set  G,  where  the  positive  orbit  of  the  system  is 
contained  in  G,  then 


if  =  <V,V>y 

if=  2<v,v>v  =  2<Av,v>v  +  2<Bu,v>v 
k  —  2<Bu,v>v  =  2<u,B*v>u 

This  leads  to  construction  of  a  feedback  control  u,  where 

u  =  — B*v 


Hence, 

if  =  -2lB*vll^  <  0 


(7-56) 


(7-57) 


(7-58) 


Since  X  is  negative  with  respect  to  the  norm  in  the  control  space  U,  no  further 
implications  can  be  made  about  the  negative  definiteness  of  if  in  terms  of  the  norms  in 
state  space  V.  This  means  that  the  Lyapunov  direct  method  is  not  applicable  to  such 
systems.  However,  from  the  invariance  principle  it  can  be  concluded  that  if  the  positive 
orbit  of  the  system  is  compact,  then  the  closed  loop  system  S(t)v  — ►  as  t  — ►  cc, 

where  M'^  is  the  largest  invariant  set  on  which  if  =  0.  The  motion  of  the  dynamical 
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system  is  asymptotically  stable  and  S(t)v  — *  0  as  t  — ►  0  if  =  {O}.  Therefore,  it  is 
sufficient  to  add  approximate  controllability  in  order  to  arrive  at  stabilizability  of  this 
system,  as  shown  by  [38]  in  the  following  theorem: 

Theorem  7.2.  Let  A  in  system  (7-24)  be  the  infinitesimal  generator  of  a  contraction  Co- 
semigroup  T(t)  on  V  for  t  >  0.  If  (i)  for  every  yeV,  S(t)y  remains  in  a  compact  set  of  V 
where  S(t)  is  the  semi-group  generated  by  closed  loop  operator  as  described  above,  and 
(ii)  (7-24)  is  approximately  controllable,  then  S(t)y  — *•  0  as  t  — ►  oc,  i.e.,  system  (7-24)  is 
weakly  stabilizable. 

Proof.  If  C  =  A  —  BB*  is  the  closed  loop  operator  with  the  control  law  (7-57)  and 
S(t)  is  the  semigroup  generated  by  C,  then  C*  =  A*  —  BB*  is  the  generator  of  S*(t)  for 
the  following  system 

y  =  C*y  (7-59) 

y(0)  =  Vo  yeV 

Consider  the  Lyapunov  functional 

^  =  <y,y>v  , 

then  — 2l[B*yllu  =  — 2l(B*S*(t)vollu,  similar  to  (7-58).  From  condition  (i)  and 

<  0,  it  can  be  concluded  that  S*(t)vo  — *•  as  t  — ♦  oc,  where  is  the  largest 

invariant  set  on  which  i?  =  0.  From  condition  (ii)  it  can  be  concluded  that  =  |0}. 
Let  meM'*'  and  define  z(t)  as 

t 

z(t)  =  /  S*(s)m  ds  (7-60) 

O 

Since  C*  is  closed,  z(t)eD(C*)  =  D(A*)  C  V,  and  z(0)  =  0,  then 
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i(t)  =  (A*  -  BB*)z(t)  +  m 


from  the  definition  of  M"'’,  it  follows  that 


(7-61) 


k  =  -2lB*S*(t)  mll?r  =  0 


(7-62) 


or 


B*S*(t)  m  =  0,  for  all  t  >  0 


(7-63) 


This  implies  that  B*[z(t)  —  z(0)|  =  0,  or 


B*z{t)  =  0 


(7-64) 


Hence, 


z  =  A*z(t)  +  m 


(7-65) 


and 


z(t)  =  /  T*(3)  m  ds 


(7-66) 


From  (7-64),  it  can  be  concluded  that 


B*z(t)  =  0  =  /  B*T*(s)  m  ds 


(7-67) 


Therefore, 


B*T*(s)  m  =  0 


-68) 


However,  if  system  (7-24),  by  assumption  (ii),  is  approximately  controllable,  then 


B*T*(s)  m  =  0  implies  m  =  0 


This  indicates  that  every  element  of  the  invariant  set  M"''  is  zero.  Therefore,  from  the 
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invariance  principle,  S*(t)V(,  — *■  0  as  t  — ♦  oo.  Thus,  S(t)vo  — *■  0  weakly  in  V  as  t  — ►  cc. 

□ 

The  analysis  in  this  chapter  indicates  that,  unlike  the  case  of  finite  dimensional 
systems  where  controllability  implies  stabilizability,  for  infinite  dimensional  systems  such 
an  implication  is  not  appropriate.  In  general  for  infinite  dimensional  systems,  exact 
controllability  provides  sufficiency  for  existence  of  a  stabilizing  control.  The  special  case 
of  an  infinite  dimensional  dynamical  system  which  has  the  characteristics  of  a 
contraction  semigroup  is  considered  in  the  last  section  of  this  chapter.  In  this  case, 
approximate  controllability  along  with  compactness  of  the  closed  loop  motion  provided 
sufficient  conditions  leading  to  the  existence  of  a  stabilizing  control.  It  can  be  shown 
that  under  the  condition  of  contraction,  the  compactness  of  closed  loop  motion  would 
imply  the  close  range  of  the  control  operator  C^.  This  in  turn  implies  that  the 
approximate  controllability  would  lead  to  an  exact  controllability. 
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CHAPTER  8  -  SUMMARY  AND  CONCLUSIONS 

A  summary  of  the  research  presented  in  the  previous  chapters  and  the  results  are 
reviewed  in  this  section.  The  major  conclusions  drawn  from  the  study  are  presented, 
and  some  recommendations  for  future  directions  are  described. 

Summary 

The  semigroup  properties  of  states  of  dynamical  systems  were  presented.  The 
extension  of  the  Lyapunov  direct  stability  method  for  distributed  parameter  systems  was 
presented.  It  was  shown  that  whenever  the  Lyapunov  functional  subject  to  the  direct 
method  or  the  equilibrium  point  does  not  exist,  then  asymptotic  behavior  of  the 
distributed  parameter  systems  can  be  predicted  by  extension  of  the  invariance  principle. 

A  simplified  model  of  the  MPD  thruster  using  conservation  laws  and  Maxwell’s 
equations  was  derived.  This  model  was  used  for  stability  and  controllability  analysis  of 
the  MPD  systems.  The  theorems  and  concepts  of  the  stability  of  DPS  was  applied  to  a 
parabolic  model  of  a  magneto-plasma  dynamic  system,  subject  to  perturbations  about  its 
null  equilibrium  velocity. 

Asymptotic  stability  of  the  transverse  and  longitudinal  modes  of  motion  were 
derived.  Similar  treatment  was  applied  to  the  case  of  the  MPD  accelerator  with  nonzero 
equilibrium  fiow  velocity.  The  stability  of  the  linearized  model  was  derived  based  on  an 
equivalent  norm.  In  addition,  the  stability  of  the  original  nonlinear  DPS  model  was 
investigated  and  derived. 

The  concepts  of  controllability  and  observability  for  linear  time  invariant  DPS 
were  compared  with  the  corresponding  concepts  for  finite  dimensional  system. 
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Fuadamental  differences  between  these  two  types  of  systems  from  controllability  and 
observability  points  of  view  were  described. 

The  stabilizability  of  DPS  in  the  absence  of  asymptotic  stability  was  analyzed  and 
the  resulting  theorems  were  applied  to  a  class  of  linear  symmetric  hyperbolic  systems.  It 
was  shown  that  if  the  system  is  exactly  controllable,  then  existence  of  a  stabilizing 
distributed  control  can  be  achieved.  For  a  special  class  of  DPS  where  states  of  the 
system  form  a  contraction  semigroup,  stabilization  was  proven  based  on  the  invariance 
principle  and  approximate  controllability. 

Conclusions 

Following  items  are  concluded  from  the  research  presented  in  this  manuscript. 

In  general,  Lyapunov  stability  theorem  provides  suflBcient  condition(s)  for  the 
stability  or  asymptotic  stability  of  systems.  In  the  case  of  systems  represented  by  a  large 
number  of  partial  differential  equations,  Lyapunov’s  method  provides  an  applicable 
process  for  stability  analysis  as  opposed  to  the  spectrum  analysis  of  high  order 
characteristic  equations  with  the  presence  of  wave  number  in  the  characteristic  equation. 
Moreover,  in  infinite  dimensional  systems,  negative  definiteness  of  supremum  of  the 
spectrum  provides  only  the  necessary  conditions  for  asymptotic  stability.  However,  the 
extra  condition  from  the  spectrum  determined  growth  assumption  should  exist  to 
guarantee  the  existence  of  asymptotic  stability. 

In  a  case  that  asymptotic  stability  does  not  exist,  then  exact  controllability 
provides  a  sufficient  condition  for  the  stabilizability  of  the  system.  I  general,  exact 
controllability  is  very  difficult  to  obtain  and  has  been  proven  only  for  a  few  special  cases. 
This  study  has  shown  that  the  stabilization  of  linear  symmetric  hyperbolic  systems  with 
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distributed  control  is  plausible  without  the  need  for  exact  controllability. 

The  results  of  this  study  can  be  applied  to  more  elaborate  models  of  MPD 
thrusters,  and  more  definite  answers  on  the  subject  of  system  observability  and 
stabilization  can  be  derived. 

Recommendations 

The  following  issues  are  related  to  the  materials  discussed  during  the  course  of  this 
research  and  are  recommended  as  topics  for  future  research  in  continuation  of  this  study 
on  distributed  parameter  systems: 

1.  determination  of  different  types  of  strategic  points  for  general  classes  of  DPS  and 
their  applications  to  sensor/actuator  optimal  locations, 

2.  closed-loop  robust  control  of  DPS  with  respect  to  noise  and  random  excitations, 

3.  effects  of  delay  in  DPS  and  robust  control  design  for  delayed  DPS, 

4.  stabilization  of  nonlinear  structural  systems  based  on  DPS  models, 

5.  experimental  investigation  of  the  proposed  stabilizability  technique, 
implementation  of  the  results  on  an  actual  MPD  thruster  at  one  of  the  Air  Force 


Laboratories 
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Appendix  A:  Proof  of  Theorem  3.7 

Given  any  e  >  0,  and  without  any  loss  of  generality,  let  e  <  r.  Defining 
a{e)  =  inf»,_Y^ii.ti?,  and  noting  that  a  >  if(v,)  +  f(c),  it  is  clear  that  the  ball  of  radius  e 
at  V,  i.e.,  Bt(v,)  contains  a  disjoint  component  of  the  set  (veV  |  i?(v)  <  a}  such  that 
V,  e  G®.  Since  i?  ;  V  — ►  R  is  continuous,  G^  is  an  open  and  positive  invariant  set  and 
there  exists  fi(c)  >  0  such  that  Bi(v,.)  C  Gq.  Therefore,  given  any  v,,  e  B5(ve)  C  Gq.,  the 
positive  orbit  is  contained  in  Ga,  i.e.,  ^(v^)  C  G^  CBj(Ve),  Therefore,  the  point  Vg  is 
stable. 

With  a  and  S  defined  as  above,  and  0  <  e  <  r,  consider  any  fixed  point 

Vo  e  Bf(Ve)  C  Gft.  Theorem  3.6  implies  that  i?(T(t)Vo)  is  non- increasing  on  R'''. 

Therefore,  as  t  — ►  oo,  i?(T(t)vo)  — ►  inf  i?(T(t)vo)  =  /^(vq)*  Clearly,  <»>/?>  ^{^e) 

UR* 

because  .2’{v)  >  i?(v,)  for  all  vcGaCBr(v,).  Either  /3  must  be  ^(ve)  or 
'y(Vo)  n  (vcGa  (  if(v)  <  0}  for  0  >  i?(ve)  is  empty.  In  the  latter  case,  continuity  of  X 
implies  that  there  exists  a  i/  >  0  such  that  7(vo)  fl  B,,(ve)  is  empty,  assuming  that 
0  ^  i?(ve)  and  i?(v)  <  — g(d(v,Ve))  for  all  veG^  C  Br(v).  From  the  hypothesis,  g(-)  is  a 
monotone  function,  hence  g{u)  >  0  and 

t 

if(ve)  <  i?(T(t)vo)  =  if(vo)  -  /  g{u) 

0 

<  if(vo)  —  t  g(i/)  — ►  — oo  as  t  — ♦  oc 

which  is  impossible.  Hence,  3  =  if(ve)  and  if(T{t)vo)  — ^(vg)  as  t  — ►  CO.  This  implies 
that  f(d(t(t)vo,Ve))  — V  0  as  t  — ►  oo.  Since  f  is  a  monotone  function,  T(t)vo  — ►  Vg  as 
t  — ►  oo  and  v,  is  asymptotically  stable  [28]. 
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Appendix  B:  Proof  of  Theorem  3.11 

If  V  is  complete  and  'y(v)  is  precompact,  then  ri(v)  is  nonempty  and  invariant; 
moreover,  d(T(t)v,fl(v))  — *•  0  as  t  — *•  oo.  Assuming  that  ^(v)  C  G,  one  has 
n(v)  C  7(v)  C  G.  If  'y(v)  is  not  precompact,  n(v)  may  hs  empty.  In  this  case,  the 
theorem  is  obviously  true  but  meaningless.  Hence,  the  case  that  r2(v)  is  nonempty  is 
considered.  Since  if{v)  <  oo  and  i?(T(t)v)  is  nonincreasing,  then  if(T(t)v)  has  a  finite 
value  for  teR'*’.  This  implies  that  if— *-/?<ooas  t— ^oo,  where  =  infteR-if(T{t)v). 
Since  fl(v)  is  assumed  to  be  nonempty,  n{v)  C  7(v)  C  G,  and  since  if(v)  is  continuous 
from  the  definition  of  n(v),  it  follows  that  ^(z)  =  /?  for  every  z  e  n(v).  Furthermore, 
since  n(v)  is  positive  invariant,  if(z)  =  0  for  every  z  e  f2(v)  and  proof  is  complete  [26]. 


Appendix  C:  Derivation  of  Model  for  MPD  Systems 


The  plasma  dynamic  equations  for  this  system  consist  of  Maxwell’s  equations, 
Ohm’s  law,  conservation  of  electric  charge,  equation  of  state  (ideal  gas  law)  and  a  set  of 
mass,  momentum  and  energy  equations  [70].  In  derivation  of  these  equations  it  is 
assumed  that  the  bulk  properties  of  the  gas  (plasma)  are  shared  by  all  species  contained 
in  the  plasma. 


Maxwell  equations: 


VX 


ae  E 
dt 


V  X  E  =  - 
V-H=0 


5/ieH 


Ohm’s  law: 


Ji  =  ctj^Ei  4- /ie(  u*  X  H)i  +  PeU*i  i  ^  direction  index  x,y, z  or  xi ,X2,X3 


Conservation  of  electric  charge: 


3Pe  ^  <5Jj  . 

— — I-  V  -r —  =0  j  ^  direction  index 
dt  (9xj 


Equation  of  state  (ideal  gas  law): 


P  =RpT 


P  =  thermodynamic  pressure 
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Conservation  of  Mass: 


Conservation  of  Momentum: 


Du*i 

Dt~ 


dPt  3  <^ii 


where, 


P  =  sum  of  thermodynic  pressure  and  radiation  pressure,  the  latter  is  negligible. 
Fe.  =  Pe  Ei  +  /ie(7  X  H)i  i  ^  X,y,Z  or  Xi,X2,X3 

Fj.  ==  gravity  force  per  volume,  negligible. 

5u*i 


Fy  =  collisional  forces,  negligible. 
Energy  Equation: 


,  3  dpe^  u*j  _  _  dn*jPi  dn*^Ti^  *  r  4-  ^ 
dt  5xj  5xj  dxj  ^  ^  dxj 

where  e^m  =  total  enprgy  per  unit  mass  cii  CyT. 

Qj  =  heat  flux  in  direction  j. 

In  this  model,  it  is  assumed  that  the  plasma  is  originally  at  rest  with  the  pressure 
Pq,  temperature  T^,  and  density  p^.  An  external  uniform  magnetic  field  Hq  is  applied  to 
the  system,  where 

Hq  =  i  Hj  +  j  Hy  +  k  0  . 
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There  is  no  electric  field  applied  to  the  system.  Plasma  is  perturbed  by  a  small 
disturbance  and,  as  a  result,  the  state  of  the  system  is  a  combination  of  the  stationary 
(equilibrium)  part  and  the  perturbed  portion.  The  velocity  vector  for  the  basic  flow  is 
zero.  Therefore, 

u*  =  i  u(x,t)  -I- jv{x,t)  +  k  w(x,t) 

However,  an  assumption  is  made  that  the  variations  of  variables  are  only  functions  of 
one  spatial  dimension,  x,  and  time.  Therefore,  instantaneous  pressure,  temperature  and 
density  can  be  written  as 

P=Po-l-P'(x,t) 

T  =  To  +T'(x,t) 

P  =  Po  +  P  (x,t) 


Electric  and  magnetic  fields  can  be  represented  as 


E  =  i  Ex(x,t)  -I- j  Ey(x,t)  +  k  E*(x,t) 
H=Ho  +h(x,t) 


Hx  +  h,(x, 


Hy  -t-  hy(x,t) 


+  k  hi(x,t) 


Current  density  J  and  net  electric  charge  Pg  are 


J  J(x,  t),  Pg  Pe(^>t)  . 

^(  •\  ^(  •\ 

The  one-dimensional  assumption  results  in:  =  0,  =  0. 

ay  oz 

Inserting  the  above  simplification  into  the  set  of  general  dynamic  equations,  the 
following  describing  equations  can  be  derived: 


Maxwell’s  equations: 
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5E, 

(C-1) 

5Ey  (9hj 

"  a  <9x 

(C-2) 

(9E-  (9hy 

“  -sr 

(C-3) 

II 

o 

(C-4) 

5hy  dEt 

;c-5) 

dht  5Ey 

dt  ~  dx 

(C-6) 

Generalized  Ohm’s  law; 

Jy  =  <7(Ex  -  Me  wHy)  +  PeU 

(C-7) 

Jy  =  <7(Ey  +  Me  WHy)  +  P^V 

(C-8) 

J,  =  a{Et  +  MeUHy  -  MeVHx)  +  PeW 

(C-9) 

Conservation  of  electric  charge: 

5pe  ^Jx 
dt  dx 

(C-10) 

The  equation  of  state  for  perturbed  variables  is 

rj'  *  rp' 

-5-  =  —  +  ,  where  =  PoRTq 

Po  Po  To 

(C-U) 

The  linearized  continuity  equation  becomes: 

dp  du 

(C-12) 

The  linearized  equations  of  momentum  are; 
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^0  $■  =  -  ^  +  4- 5- - J- +  "•E. 


dt 


dv 


dx 

c^v 


a  ^  dx^ 


3  dx^ 
f^e  Pe  ®y 


M  -^  +  /^(Ji  H,  -  H,  Jy)  +  P.E. 


(C-13) 

(C-14) 

(C-15) 


It  is  assumed  that  the  nonliiiear  perturbation  terms  are  negligible  in  comparison  with 
the  linear  terms.  Therefore,  the  energy  equation  becomes 


PoC 


^  =-p^RT„-^+K 


OV-v  ^  ^ 


d^T' 

dx^ 


(C-16) 


Decoupled  Modes  of  Motion 


In  the  case  of  a  neutral  plasma,  i.e.,  Pe  —  number  of  ions  and  electrons  per 

volume  of  plasma  are  nearly  equal.  For  this  case,  if  one  considers  the  fact  that 
5hx  (9hx 

—  0,  =  0,  then  it  is  possible  to  distinguish  between  two  modes  of  wave 

ot  ax 

propagation:  transverse  mode  (z-direction)  and  longitudinal  mode.  In  the  transverse 
mode,  the  states  are  found  to  be  h,  and  w,  and  the  state  equations  can  be  formed  from 
the  reduction  of  equations  (C-l),  (C-2),  {C-6),  (C-7),  (C-8),  (C-10),  and  (C-15).  The  rest 
of  the  equations  can  be  reduced  to  the  form  state  equations  for  the  longitudinal  mode. 


(i)  Transverse  Mode 

The  magneto-gas-'^ynamic  assumption  results  in  an  insignificant  magnetic 
induction  effect  in  Maxwell’s  equations  from  the  terms  carrying  variations  of 
electric  field  with  time.  This  is  due  to  the  fact  that  nondimensional  parameters 

E„ 


toU* 

~e~ 


and  Re  = 


are  of  the  order  of  one  or  smaller,  and 
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Rc  =s  —5-  =  <K  1  [70].  The  resulting  equations  are 

c* 

Sbf  ^h,  ^ 

"ar  " ar 

<9w  _  d^w  V*  (9hi 

~  hT  ^ 


(C-17) 

(C-18) 


where  — -  and  V, 


The  parameter  Vj^  is  defined  as  the  x- 


component  of  the  speed  of  the  Alfven  wave. 


(ii)  Longitudinal  Mode 


The  state  equation  for  this  mode  can  be  reduced  to 


<9hy 


=  i/]g- 


H  -l-H 

-H, 


5x* 


dy  _  <9*v  Vx  5hy 
di  dy^  Hx  ^ 


-  RTc 

a  ' 


■  ar" 

dx  dx 


4  5*u 


V2y 


(9h„ 


Hv  a 


a  a 

ar"  ^  K  ^t"  _  ^  ^ 

PoC./  dx^  Cy  ^ 


(C-19) 

{C-20) 

(C-21) 

(C-22) 

(C-23) 


where 


n 

p 


The  parameter  Vy  is  defined  as  the  y- 


component  of  the  speed  of  the  Alfven  wave. 


